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SYNTHESIS PROCEDURE FOR LINEAR TIME-VARYING 
FEEDBACK SYSTEMS WITH LARGE PARAMETER IGNORANCE 


Abstract- -This work is addressed to the development of synthesi 
procedures for linear time-varying feedback systems. It is 
assumed that the plant can be described by linear differential 
equations with time-varying coefficients; however, ignorance is 
associated with the plant in that only the range of the time- 
variations are known instead of exact functional relationships. 
As a result of this plant ignorance the use of time-varying 
compensation is ineffective so that only time-invariant com- 
pensation is employed. In addition, there is a noise source 
at the plant output which feeds noise through the feedback 
elements to the plant input. Because of this noise source the 
gain of the feedback elements must be as small as possible. 

No attempt is made to develop a stability criterion for time- 
varying systems in this work. 

Two synthesis procedures are developed and investigated. 
Both procedures assume system specifications to be given in 
the frequency domain (i.e., specifications on the Fourier trans 
form of the system signals) and both procedures arrive at the 
transfer function for the linear compensation. The method used 
in the procedures takes an iterative approach in that a design 
is made, then the system is simulated on a computer to test its 
adequacy. If the design is unsatisfactory, a subsequent design 
is made until either the system is satisfactory or instability 
is reached. Application of the procedures to general systems 
is limited because of the lack of a stability criterion; how- 
ever, the procedures can be successfully applied to systems 
with time-varying gains since a stability criterion exists 
for such systems . 
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CHAPTER I 


INTRODUCTION 

1.1 Description of the System 

A block diagram of the system with which this work is con- 
cerned is shown in Figure 1.1. 



Figure 1.1 

Block Diagram of the Feedback 
System Under Study 

i 

The block P represents the plant which can be described by a 
linear differential equation with time-varying coefficients and 
blocks G and H represent compensations. There is also associated 
with the plant a form of ignorance in that it is assumed the 
time— varying coefficients have known bounds on their variations 
but the exact functional relationships of the variations are 
unknown. All components are considered to be single input 
single output. The compensations are to be designed so that the 
desired system response is obtained while the effects of the 
plant time variations on the system output are reduced to an 
acceptable level. This is to be done with a minimum amount of 
noise transmitted to the plant input due to N. 
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1.2 Histroical Background 

Prior to 1950 there appeared to be little interest in time- 
varying control systems, and consequently little research had 
been done on such systems. * However, since 1950 more interest 
has been shown in time-varying control systems primarily due to 
the advent of such modem technology as rocket flight, space 
exploration, and control of high performance aircraft. The 
behavior of the plants of these systems can often be approximately 
described by linear differential equations with time-varying 
coefficients; thus, techniques have been sought which will aid 
the control engineer in designing a suitable system around a 
linear time-varying plant. ^ 

The major portion of the work to date on linear time- 
varying control systems has been done either directly in the 
time domain, or some type of transform has been defined with the 
subsequent investigation of this transform. 

Time Domain Approach : 

The time domain approach is the study of the system behavior 
as functions of time. It is well-known that the output of a 
linear system y(t) is related to the input of the system x(t) by 
the superposition integral 

y(t) = / fc <j>(t,x) x(t) dr ( (1-1) 

—CO 

where <j>(t,x) is the output of the system due to a unit impulse 
applied at t = t and is defined as the system impulse response. 

For the time invariant case <j>(t,x) is a function of t - x. 



Friedland^ has formulated a representation of the impulse response 
<f>(t,r) as a lower diagonal matrix where the i, jth element in the 
matrix is the value of <j)(t,T) at time t = t^ due to an impulse 
applied at time t = t j . This type of representation lends itself 
readily to the analysis of sampled data systems since at the 
sampling instances the integration is replaced by matrix multi- 
plication. Cruz^ has worked bn the synthesis of time-varying 
control systems using the Friedland representation. The compen- 
sations which he obtains are in the form of impulse responses. 
These compensations are inevitably time-varying themselves, which 
leads to the study of the synthesis of prescribed impulse 
responses of time-varying systems . 1 

An algebra has been presented independently by Darlington,® 

A. V. Solodav,^ and Stubberud^O which allows block diagram 
manipulations of a system represented by differential equations. 
This would allow the blocks of the system shown in Figure 1.1 
to be combined into a single block represented by a linear time- 
varying differential equation and could then be more easily 
analyzed. Stubberud also recognized that if the system specifi- 
cations were given in the form of a desired differential equation 
of the overall system, then the method could be used as a synthe-t 
sis procedure. He has thus presented two synthesis techniques 
for synthesizing a desired system differential equation from the 
time-varying differential equation describing the plant. As in 
the techniques employed by Cruz, the compensations which are 
obtained by Stubberud are inevitably time-varying. 
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Transform Approach: 

Largely due to the success with which the Laplace and Fourier 
transforms can be applied to differential equations with time 
invariant coefficients, a number of investigators have sought 
transform techniques which can be applied with equal success to 
linear differential equations with time-varying coefficients. 

The desired goal in these investigations is to arrive at a 
technique which will yield general solutions to linear time- 
varying differential equations merely by solving an algebraic 
equation in the transform domain in much the same manner as a 
differential equation with constant coefficients can be solved by 
using the Laplace transform. Unfortunately, up to the present 
time no transform technique has been found that can be applied to 
a linear time-varying system in a practical manner. 

Transforms have been found for special types of time-varying 
differential equations, These special integral transforms include 
the Laplace transform, the Mellin transform, and the Hankel 
transform . 11 ’ 12 The Laplace transform is used in the solution of 
linear differential equations with constant coefficients, while 
the Mellin transform is used to solve the equidimensional 
(Euler-Cauchy) differential equation which is of the form 

t n d n y(t-) + a t n_1 ..^ n . L-y(t) . . . + a Q y(t) = f(t), 

dt n n_1 dt n_1 

and the Hankel transform can profitably be applied to Bessel's 
equation which is of the form 



t 2 dZ y.(ti + t Ml> + (at 2 - n 2 ) = f (t) . 
dt 2 dt 

Unfortunately, neither the equidimensional equation nor Bessel's 
equation possesses sufficient generality to be applied to the 
time-varying system synthesis problem. 

Aseltine 1 has proposed a method of deriving a compatible 
transform for a linear time-varying differential equation; 
however, practically, the transform can only be obtained for a 
second order time-varying system. Further, a block diagram 
algebra cannot be developed since a transform which is compatible 
with one part of a system will not generally be compatible with 
another part of the system. 

Naylor 1 ^’ 15 has taken an interesting approach in attempting 
to obtain a system transformation. Working with the Friedland 
characterization, he has defined the system transformation as a 
matrix which, when properly applied to the Friedland characteriza- 
tion of a system, results in a system matrix which is diagonal 
and is thus easily manipulated. However, the transform domain of 
the Naylor Transformation does not have a simple interpretation, 
so that synthesis in this transform domain appears to be as 
difficult as working directly in the time domain. 

L. A Zadeh 16,17,18,19 has defined a system transformation 
that is similar to the transfer function of a time-invariant 
system. The transformation proposed by Zadeh which is referred to 
as the system function is defined as 
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<f(t,ja)) =» / tj>(t,x)e T ^dx 

""00 

where <f) ( t , t ) is the impulse response of the system. Cruz^O has 
made a study of the system function and has proposed a technique 
for synthesizing a desired system function provided, of course, 
that such a desired system function is known. This work does not 
consider the feedback problem. Kaplin^ has considered the system 
function from a more analytical viewpoint and has presented 
theorems which give some of its mathematical properties . Other 
researchers^ ,23 ,24 , 25, have also considered the Zadeh trans- 
formation, but no significant progress has been made past the 
purely mathematic description and the study of its general 
characteristics. 

I 

1.3 General Considerations 

In this section three of the more important aspects which will 
affect the design of the system described in Section 1.1 will be 
considered. 

Plant Ignorance : 

Although the plant can be described by a linear differential 
equation with time-varying coefficients, it is assumed that these 
variations are not known explicitly as functions of time. That 
is, the plant has associated with it some type of ignorance. Such 
an assumption is necessary in the formulation a'f the feedback 
problem; otherwise, it would be possible in theory, at least, to 
design a prefilter which would effectively cancel the parameter 



variations of the plant and feedback would not be mandatory 26 . 

It will, therefore, be assumed that ranges for the values and the 
rates of variations of the time-varying coefficients are known 
but that the time of occurrance of the variation is unknown and 
the exact rate and values of the variations are unknown. Figure 
1.2 gives an illustration of what is assumed to be known and 
unknown about the coefficients. Referring to Figure 1.2, the 
range of a n (t) is assumed to be A < a n (t) < B, and the range of K 
is assumed to be Rj < K < K 2 ; however, the actual values of A', 
B', and K as well as x are assumed to be unknown. 



Figure 1.2 

Illustration of the Ignorance Assumed to 
Be Associated with Each Coefficient 

Due to the assumption that the plant is not completely known 
it is evident that the synthesis procedures of Stubberud 10 and 
Cruz, 7 which arrive at time-varying compensations as mentioned in 
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! Section 1.1, cannot be used since these procedures rely on knowing 

the equation coefficients explicitly as functions of time. Never- 
theless, one may still consider whether or not time-varying 
compensation might be superior to time-invariant compensation. 

Under the ignorance assumptions just outlined, time-varying 
compensation has no particular advantage over time-invariant 
compensation. The argument leading to this conclusion is as 
follows. The compensation cannot be used to cancel out 
the plant parameter variations since they are not known exactly; 
therefore, if time-varying compensation has an advantage over 
time— invariant compensation, it must be due to the time— varying 
compensator's superior signal processing properties or filtering 
properties. However, to take advantage of the time-varying 
compensator's signal processing properties, the time-varying 
characteristics of the signal to be processed must be known. 

But, these are not known in general due to the assumption of 
plant ignorance which thus leads to the conclusion that there is 
no particular advantage in using time- varying compensation. A 
second consideration from a more practical point of view is the 
fact that a time-varying compensator is much more difficult and 
expensive to fabricate and wauld in all probability be physically 
larger than a time-invariant compensator. Thus, even when the 
plant is known exactly, it may be more advantageous to use time- 
invariant compensation. For these reasons only time-invariant 



compensations will be considered in the synthesis of the system 
of Figure 1.1 

Specifications : 

Several performance criteria have been established for linear 
time- invariant systems which include step response criteria, 
steady state error criteria, stability criteria and frequency 
response criteria. ^7 Because of the linearity of the system, a 
modified form of the majority of these criteria can still be used 
in the linear time-varying system. Unfortunately, the stability 
criteria of a linear time-invariant system cannot be applied to 
a linear time-varying system. 

The design criteria which will be used in this work is the 
step response of the system. To justify the use of the step 
response to judge system performance of linear time-varying 
systems, one can argue that any input signal can be approximated 
by a sum of step functions; and, due to the linearity of the 
system, the system output will be the sum of the step responses 
due to each individual step function in the input. Thus, a 
reasonably well-behaved step response will imply a reasonably 
well-behaved response to any input. A second justification 
can be made by examining the reason for using the step response 
as a performance criteria in a linear time-invariant system. In 
a linear time-invariant system the impulse response <J> (t - x) and 
the step response c(t) are related by 
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4>(t - T) = i_ C(0 

dC 


5 = t - T 


and knowing <|>(t - t) the system response to any input can be 
determined. Thus, the performance of a linear time-invariant 
system can be determined from its step response. If one could 
obtain the impulse response of a time-varying system from its 
step responses, then it would be possible to associate system 
behavior with these step responses in the same way as is done in 
the time-invariant system. Note that in time-varying systems 
one would actually consider a number of different step responses, 
each step response corresponding to a particular time of interest 
It will now be shown that the impulse response for a linear 
time-varying system can be constructed from a set of step 
responses of the system. Suppose the values of all system step 
responses starting in the interval from t = 0 to t = tj. are 
known at t = t 1 . Figure 1.3 shows one such step response c y (t) 
which is due to a step input applied at t = y. The value of a 
step response evaluated at t ■ tj due to a step input applied at 
t = y is given by 

/tl 

c y (t 1 ) = J (j) ( t x , x) dt . 


( 1 - 2 ) 
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Figure 1.3 

Step Response Due to a 
Step Input Applied at t = y 


By considering t^ fixed and varying y, c^(t^) becomes a function 
of y. Let 


A(y) = - c y ( tl ) 

so that Equation (1—2) can be written 

f Y 

A(y) = J t <f> (t 1 ,T)dx. 

Differentiating X(y) with respect to y one obtains the desired 
result 


dX(y) _ 

dy 


<t> (t l5 y) . 


Thus, with a time-varying system it is possible to associate 
system behavior with its step responses. 
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The design specifications will be assumed to be given in the 
form of an acceptable step response. In principle step response 
specifications can be given either in the time domain or the 
frequency domain. Generally , step response specifications are 
given in the time domain since one is normally aware of the 
desired time domain response. The difficulty with frequency 
domain specifications lies in the fact that it is not known 
precisely what constitutes acceptable specifications. The basic 
problem of translating time domain specifications into frequency 
domain specifications and vice versa has yet to be solved^ , 
although investigations of single responses have been made by 
researchers in the field . 29,30 Nevertheless, frequency domain 
specifications are desirable since by using such specifications 
general design procedures can be developed which are not dependent 
upon the order of the differential equation describing the system, 
and thus such procedures are not limited to simple systems. 

Further, the "cost of feedback" is clearly evident in the fre- 
28 

quency domain . 

This study will not attempt to solve the problem of translat- 
ing time domain specifications into frequency domain specifica- 
tions; rather, the necessary frequency domain specifications will 
be assumed known. These specifications will be of the form of 
a limit on the magnitude of the Fourier transform of the 
difference in the actual system step response c(t) and a desired 
system step response c Q (t). Some justification for such 
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specifications can be seen from the Fourier transform of this 
difference 

E(joi) = / e(t)e - 3 wt dt 
—00 

where 

e (t) = c ( t) - c Q (t) . 

If E(jcu) is "small" over the entire frequency range then e(t) 
will also be "small". The exact relation between the magnitude 
of the argument of E(jm) and the magnitude of e(t) has yet to be 
determined . 

Stability Considerations: 

Although a number of researchers have investigated the 
stability of linear time-varying systems, necessary and sufficient 
conditions have yet to be found which will insure the stability 
of the system of Figure 1.1?* A number of definitions and con- 
cepts of stability exists; however, in this work a system will be 
considered stable when it is bounded-input bounded-output stable. 
That is, if the system input is bounded, then the system output 
will also be bounded for a stable system. 

Research in the stability of a time-varying system has taken 

basically two approaches. The first is to examine the system's 

32 33 

differential equation or impulse response in the time domain; * 
however, no practical results have been obtained which lend 
themselves directly to the application to the system of Figure 1.1. 
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The second approach has been to obtain a frequency domain sta- 
bility criterion similar to the Nyquist stability criterion of 
time-invariant systems 5 33 Results have been limited to systems 
in which the only time- varying element is the plant gain. Such a 
system is illustrated in Figure 1.4. 



Figure 1.4 

System with Time- varying Gain 

However, even for this case only sufficient conditions for 
stability have been determined rather than the necessary and 
sufficient conditions which would be required for a completely 
satisfactory stability criterion. 

This thesis will not attempt to develop a stability criterion 
for linear time-varying systems; rather, synthesis procedures will 
be developed under the assumption of a stable system. It will be 
seen that the procedures can easily be modified to include any 
applicable stability criterion that specifies constraints on the 
transfer function of the feedback compensation H. The sufficient 
conditions for stability which have been developed for the system 






of Figure 1.4 represent such a criterion and the synthesis pro- 
cedures will be modified appropriately to include this criterion 
in the design of systems of the type shown in Figure 1.4. 



CHAPTER II 


DEVELOPMENT OF INPUT-OUTPUT RELATION AND SYSTEMS EQUATIONS 
2.1 General 

A number of representations exist which relate the output of 
a linear system to its input. The differential equation itself 
can be considered as one such representation. Unfortunately, the 
differential equation is not directly useful as a design tool 
because it does not allow for the system characteristics to be 
readily observed in order that a suitable design can be developed. 
It is thus desired that an input-output relation be found which 
is suitable for the development of a design procedure. 

An investigation of input-output relations was made in an 
effort to find one which is well suited for design purposes. The 
investigation centered on four functions in particular. The four 
functions which were considered are: 

(1) The impulse response cf> ( 1 9 x ) which is defined as the 
output of the system at time t due to a unit impulse 
applied at time x.^l 

(2) The system function which was introduced by Zadeh^ 
and is defined as 

$(t,jw) = / <f>(t ,x ) e T ^dx 

—CO 

(3) The complementary function which is defined by 

$(jm,x) = / <j>(t s x)e ^ t-T ^dt 

T 
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(4) The frequency impulse response which is defined by 

OO t 

$(ju,jy) = / / <f> ( t , x )— ^ ■■ - e~^ t0t dx dt 

— oo —oo 2 7T 

is the Fourier transform of the output at frequency jco 
due to the input of an unit impulse of frequency at _jy. 
A good deal of time was devoted to the study of these four 
functions and a number of interesting properties were found and 
investigated. However, none of the functions were found to lend 
themselves to the development of a. practical design procedure. 

The difficulties encountered ranged from not being able to 
evaluate the functions to not being able to determine the output 
in terms of the input for systems having a feedback loop. An in- 
depth discussion of this study will not be presented because the 
results which were obtained generally have little bearing on the 
design problem. 

The frequency impulse response was found to be useful because 
of its notational convenience when studying the Fourier transform 
of the various signals of the system. Because it is used later in 
the work, the general characteristics of the frequency impulse 
response will be presented and a discussion of the difficulties 
in employing the frequency impulse response will be given. The 
difficulties encountered with the frequency impulse response are 
representative of the type of difficulties encountered with the 
other three. 

A clarification of notation should be made before going 
further. The term "frequency domain" refers to expressions as 
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functions of frequency such as the Fourier transform representa*- 
tion of signals, the transfer function of a system or the frequency 
impulse response of a system. The term "time domain" refers to 
expressions as functions of time such as the impulse response or 
a signal expressed as a function of time. 

2.2 Properties of the Frequency Impulse Response 

It would be possible to relate the frequency domain represen- 
tation of the system output to that of the system input by first 
employing a time domain relation such as the impulse response and 
then taking the necessary Fourier transforms. However, it is more 
convenient to employ the frequency impulse response since it is a 
direct relation between the frequency domain representation of 
the system output to that of the input. The frequency impulse 
response is defined as 

00 t . g 

‘KjWjjy) = J :f <t>(t,T) — e”^ wt dx dt (2-1) 

_oo —oo 2 7T 

where c|>(t,x) is the system impulse response. This function was 
first mentioned by Zadeh who referred to it as the bi-frequency 
transform;^ however, Zadeh gave no results of any study he may 
have made of this function. Kailath^ uses similar forms in his 
study of time-varying filters. The author feels justified in 
referring to this function as the frequency impulse response since 
it possesses many characteristics in the frequency domain that the 
impulse response exhibits in the time domain. 
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Figure 2.1 is an illustration of a linear time-varying 
system having an impulse response denoted by <f)(t,T) and a fre- 
quency impulse response denoted by $(jw,jy). The differential 
equation describing the system is assumed to be of the form 

a n (t) — x (t) + . . . + a 0 (t) x (t) = 
dt n 

b m (t) ^_y(t) + . . . + b D (t)y(t) (2-2) 

dt m 

where a n (t) ^ 0 and n > m. 


) 

Figure 2.1 

Block Diagram of a Linear System 

The primary reason for referring to $(jw,jy) as the frequency 

impulse response is because it is the Fourier transform of the 

system output evaluated at frequency jco due to an impulse of 

frequency at jy. This, of course, is the frequency domain 

equivalent of the impulse response. The property can be easily 

demonstrated as follows. Let the system input be given by 

y(t) = e ^ t which is an unit impulse at frequency jy in the 
2tt 

frequency domain^ 9 The system output is given by 3 ® 
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t -j y-j- 

x(t) = / 4><t 9 t) dT. 

—c» 2 7T 


Taking the Fourier transform of x(t)» one obtains 

«° t 

X(jw) = / / <|3 (t 9 t> — HI - dx e - ^ 1 " dt 

■ -CO -00 2 IT 


which is the definition of the frequency impulse response given in 
Equation (2-1). 

If the system is time-invariant , the frequency impulse 
response is given by 

< Kjm,jy) = $(jw) 6 (j y - jm) (2-3) 


where <£(jw) is the transfer function of the system. This result 
is obtained as follows. Since the system is time- invariant , the 
system impulse response is given by $(t s x) = <j>(t-x). Letting 
x = t-u the frequency impulse response can be written as 


$(jm,jy) 


1_ 

2ir 


/ 


- CO 

/ ( u) e~j YU du 

- —00 


e“( jw-jy) t dt. 


The integral inside the brackets is the system transfer function 
so that one can write 

< Kjw,jy) = — / <J>(jy) e d t 

2t r loo 


or upon performing the integration^ one obtains 


KjWjjy) = 0(jy) S(jio - jy) . 
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Another important property of the frequency impulse response 
is that the Fourier transform of the system output X(jw) is given 
by the relation 

00 

X(jm) = / $(jw, jy) Y(jy) dy. (2-4) 

— CO 

This expression is obtained by taking the Fourier transform of 
the output as 

00 

X(jm) = I x(t) e - ^*- dt. (2-5) 


Expressing x(t) by 


x(t) 


t 


= / 4>(t , t) y(x) 


dr 


and expressing y(x) as the inverse Fourier transform of Y(jto), 
Equation (2-5) can be written 


X(jm) 



(t > x) 


1_ 

2tt 


00 

/ Y(jy) e~j^ T dy 
—00 


dx e“j a)t dt . 


Interchanging the order of integration, the relation becomes 


X(jco) = / 


« t 


1 —jf <j>(t,x) e~ IYT e - ^ 631 " dx dt 
2 IT — oo —oo 


Y(jy) dy. 


Recognizing the term inside the brackets as the definition of the 
frequency impulse response establishes Equation (2-4) . The 
frequency impulse response can be thought of as an operator which 
operates through the relation expressed in Equation (2-4). Thus, 
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Equation (2-4) can be referred to as $(jto,jy) operating on Y(jy) 
with the result of the operation being X(jto). 

Special consideration should be given to Equation (2-4) 
because this relation demonstrates a basic difference in time- 
varying systems and time-invariant systems. Suppose the system 
were time-invariant with transfer function $(jw). Equation (2-4) 
becomes the well-known relation 

• X(joi) = 3>(ju>) Y(joi) . 

Observe that in the time-invariant system the value of the output 
at frequency oj is due only to the value of the input at m which 
is weighted by the value of the transfer function at frequency u). 
Thus, to determine the value of the output at a given frequency, 
one only needs to determine the value of the input and the 
transfer function at that frequency. For the time-varying 
system on the other hand, Equation (2-4) shows that the system 
output at a) is not only influenced by the system input at w. but 
also by all other frequencies in the system input. The amount of 
influence each frequency has on the system output is determined 
by the weighting of the frequency impulse response at each 
frequency. Therefore, to evaluate the output of a time-varying 
system, one must determine the values of the system input at all 
frequencies and the corresponding values of the frequency impulse 
response instead of just the two values required in the time- 
invariant system. The basic difficulty in time-varying systems is 
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that a much larger number of frequencies must be considered than 

in the time-r invariant case. 

If two systems are considered in parallel as shown in 

Figure 2.2, the frequency impulse response for the parallel 
combination is the sum of the frequency impulse responses of the 
two systems. This is easily obtained from the fact that the 
impulse response of two parallel systems is the sum of the impulse 
responses of the two parallel systems. 






■W — 

• ^ 

$ e = ^1 + $2 



Figure 2.2 

Combining Systems in Parallel 

Thus, the equivalent impulse response of the parallel system 
shown in Figure 2.2 is given by 

<|> e ( t , x) = ♦ 1 (t,T) + <l> 2 (t,T) 
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Applying the definition of the frequency impulse response, 

Equation (2-1) , to this relation yields the same property for the 
frequency impulse response . 

Cascade systems, unfortunately, do not combine as easily as 
do parallel systems. The equivalent impulse response of two 
systems in cascade as shown in Figure 2.3 is given by 

00 

$ e (j“>jY) - / $ 2 (jw,jD Mj5,jy) ds (2-6) 





Figure 2,3 

Combining Systems in Cascade 

This relation can be obtained by letting Y(jA) = S(jA-jy). Since 
the input is a unit impulse in the frequency domain, the Fourier 
transform of the system output is the frequency impulse response 
of the system. From Equation (2-4) the Fourier transform of the 
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system output is obtained from the intermediate signal Z(jco) by 
the relation 

CO 

X(jw) = $ e (ju),j Y ) = f $ 2 Cjto,jO Z(K) ds (2-7) 

—00 

and Z(j£) is given by 

CO 

z(j£) = / ^(jS.jD Y(jA) dA. 

—00 

Substituting Y(jX) = <5(jX-jy) and performing the integration, one 
obtains 

Z(jO = ^(jS, jy) • 

Substituting the above relation into Equation (2-7) gives the 
desired expression, Equation (2-6) . Again, thinking of the 
frequency impulse response as an operator, the cascading of two 
operators results in a new operator $ e which is obtained by 
letting the second operator $2 operate on the first operator 
Note that $ 2 operating on is not the same in general as 
operating on $ 2 . 

Cascading a time-invariant system with a time-varying system 
is quite simple. First assume that $ i n Figure 2.3 is time- 
invariant. From Equation (2-3) the frequency impulse response of 

$ 2 is 

$ 2 = $2 0“) <S(jy-jaj). 
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Substituting this expression for $2 into Equation (2-6) , the 
equivalent frequency impulse response is found to be 

$ (jw.jy) = 

Likewise, if 3^ is time-invariant, is given by 

Thus the time-invariant portions of a time-varying system are 
easily handled by ordinary multiplication rather than the more 
difficult integration process. 

9 

The concept of the inverse system is also important in 
using the frequency impulse response. Figure 2.4 shows an 
illustration of a system represented by <p and its inverse 
represented by <f> * . 



= Y 


System 


Inverse 

System 


Figure 2.4 

Illustration of System and Inverse 

The cascading of a system with its inverse results in a system 
in which the output is equal to the input. Such a system will be 
referred to as the unit system. The frequency impulse response 
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of the unit system is found from the definition to be 

= <5(jw - jy). 

Thus, if either the inverse operates on the system or if the 
system operates on the inverse, the resulting frequency impulse 
response is <5(ju) - jy). 

2.3 Discussion of the Frequency Impulse Response 

In the early stages of the investigation it was hoped that 
the frequency impulse response could be represented explicitly as 
a function of frequency and could thus then be utilized in much 
the same manner as is the transfer function for time- invariant 
systems. This hope was not realized for all efforts to obtain an 
explicit functional representation were, in general, unsuccessful. 
However, in some instances it is possible to represent the 
frequency impulse response by an explicit function of frequency. 

Any attempt to obtain the frequency impulse response from a 
direct application on the defintion is impossible except for 
systems which can be described by first order differential 
equations or by simple time-varying gains. This is due to the 
fact that the impulse response of higher order time-varying 
systems cannot be found in general. Also, for systems which can 
be described by first order differential equations a closed 
functional form cannot be obtained even with the time-varying 
function specially chosen to aid the calculation. 



28 


The simplest time-varying system consists of a time-varying 
gain f(t). The frequency impulse response for this system is 
obtained from a direct application of the definition, Equation 
(2-1) , and is found to be 

2 IT 


where F(jA) is the Fourier transform of f(t). 

The next type of system to consider is one which can be 
described by a first order differential equation of the form 

x(t) + f(t)x(t) = y (t) (2-8) 


where y(t) is the system input and x(t) is the output. The 
impulse response for this system is given by 


4> (t ,t) 


exp 


t 

- / f(t) dt . 

T 


The frequency impulse response is found by applying Equation (2-1); 
however, it is evident that unless f(t) is selected very carefully 
the integrals cannot be evaluated and an explicit functional 
representation of the frequency impulse response is not possible. 

As an example, consider a system which is described by the 
first order differential equation given in Equation (2-8) where ' 


f(t) 



+ e 
+ e 


t 

t 0 


The graph of f(t) is shown in Figure 2.5. 
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f(t) 



Figure 2.5 

Graph of f(t) Used in Example 

It is easily verified that the system impulse response is given by 

4> (t , T > = - 27 - + - U(t - t) . 
e + e 

A straightforward application of Equation (2-1) yields for the 
frequency impulse response 

$(ju),jy) = ~ J o y) + / — e -(jw-jy)t dt 

3Y + 2 (jY+1) (jaH-2) e 2t +e t 

Note that even though the function f(t) was especially chosen, a 
closed functional form of the frequency impulse response is still 
not possible. 

If the function is described by a differential equation of 
order higher than one as given in Equation (2-2), the impulse 
response cannot be determined in general and thus an expression 
for the frequency impulse response cannot be found. 



It is possible, however, to determine the inverse frequency 
impulse response of a differential equation of the form 

a n (t)~ x(t) + . . . + a (t)x(t) = y (t) (2-9) 

where y(t) is the system input and x(t) is the system output. 

The inverse system is also described by Equation (2-9); however, 
in the inverse system x(t) is the input and y(t) is the output. 
The frequency impulse response of the inverse system can be 
found by employing the property that it is the Fourier transform 
of the system output evaluated at jm due to. an impulse of fre- 
quency at j y applied to the input. An impulse of frequency is 
applied to the input by setting 


x(t) = 


_ e 3Yt 


2tt 


The output is easily found to be 

y(t) = |"a n (t)(j Y ) n + . . . + a (t)l 

*• 0 J 2tt 

Taking the Fourier transform of y(t) results in 
$ =T[y(t)J = / y (t) e" jtot dt 


jyt 


( 2 - 10 ) 


1_ 

2tt 


^(jw-jy) (jy) n + . . . + A (jw-jy) 


( 2 - 11 ) 


where A ± (jA) is the Fourier transform of a i (t) . Presumably the 
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frequency impulse response can now be obtained by finding the 
function $(jw,jy) which satisfies the equation 

00 

<$(jw - jy) = / $(jw,jx) $ _1 (jx,jy) dx. 

—00 

Unfortunately, it was not possible to find the proper function 
$(j“>jX) which satisfies the above relation for a general 
4> _1 C j x , jy) . 

Note that the right hand side of Equation (2-9) can contain 
a first order derivative of y(t) and the frequency impulse 
response of the inverse system can still be found since the first 
order differential equation can be solved. Thus, an expression 
for the inverse frequency impulse response of a system of the 
general form shown in Figure 2.6 can be obtained although the 
expression will generally contain integrals which cannot be 
evaluated in a closed functional form. 



Block 1 Block 2 


Figure 2.6 

Type of Plant for Which an Expression of the Inverse 
Frequency Impulse Response Can Be Obtained 
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It should be emphasized that if the differential equation in 
block 1 of Figure 2.6 is of order higher than one, an expression 
for the inverse frequency impulse response cannot be obtained. 

To summarize, it is seen that an expression for the frequency 
impulse response cannot be obtained for all types of systems. 
However, an expression can be obtained for time- varying gains and 
for systems that can be described by first order differential 
equations. Also, an expression for the inverse frequency impulse 
response can be obtained for a system of the type shown in 
Figure 2.6. 

2.4 Development of System Equations 

The system equations will now be developed utilizing the 
frequency impulse response. The desired objective is to develop 
an expression which will allow one to determine the Fourier 
transform of the system output from the Fourier transform of the 
system input. Equations will be found which give the desired 
relation, but, unfortunately, cannot be solved to yield the values 
of the system output. 

For convenience the system is shown in Figure 2.7. A 
detailed discussion of the noise transmission from the noise 
source N is given in Section 4.7 and will not be considered at 
this point. It is assumed that the time- invariant prefilter has 
a transfer function G(jco), and the time-invariant feedback element 
has a transfer function H(ju>), while the time-varying plant has a 
frequency impulse response P(jco,jy). 
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Feedback 

Figure 2.7 
System under Study 


One form of the basic equation is developed as follows. The 
system output is given in terms of the plant input by the 
expression 

CO 

c(jw) = / P(ju),jy) Z(jy) dy ( 2 - 12 ) 

-OO 1 

and the plant input is given by the expression 

Z(jy) = G(jy) R(jy) - H(jy) C(jy). 

Substituting Z(jy) into Equation (2-12) one obtains 
00 00 

C(jio) = / P(jm,jY)G(jy)R(jy)dy - / P( jto, jY)H(jY)C(jy)dy. (2-13) 

—00 — oo 

This is a potentially useful relation since it relates the system 
output to the system input. However, the equation is known as a 
singular integral equation and, unfortunately, no general methods' 
of solution are known for such equations. ^ Obtaining a solution 
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is doubly difficult in this case because of the complex expres- 
sion for the frequency impulse response. 

Equation (2-13) can be put in an alternate form by operating 
on both sides with the frequency impulse response of the inverse 
plant The resulting equation is 

CO 

/ P _1 0,jy) C(jy) dy = G(joj) R(jw) - H(ju>) C(jw) (2-14) 

—CO 


Again, this equation is a singular integral equation and cannot in 
general be solved. However, because of the less complex expres- 
sion for the inverse frequency impulse response of the class of 
plants which can be described by Equation (2-9), one can come 
closer to obtaining a solution from Equation (2-14) than from 
Equation (2-13) . As an example, assume the plant is a time- 
varying gain f(t), the inverse of which can be described by a 
function of the form 

f ^(t) = K + Ae“bt u(t) . 


Applying the results obtained in the last section, the inverse 
frequency impulse response is found to be 


P 1 (jw,3Y) 


K. fi(ju-jy) + — 
2tt 


A 

(b + jio - jy) 


Substituting this expression into the integral term of Equation 
(2-14) , one obtains 
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00 CO 

/ P _1 (jw,jY) C(jy) dy = J A 6(joj-jy) C(jy) dy 

— OO — OO 


00 



A 

(b-jw-jy) 


C(jy) 


dy 


The first integral on the right hand side is equal to C(jco). 

With a proper interpretation the second integral can be evaluated 
by contour integration. Instead of integrating the real variable 
from minus infinity to plus infinity, the integral is first 
multiplied and divided by j , then the complex variable s is 
substituted for the variable jy with the path of integration along 
the imaginary axis. Since C(s) is the Laplace transform of the 
system output, | C (s ) | goes to zero as |s| goes to infinity, and 
the integral can thus be thought of as being integrated around a 
closed path where the path encloses the right half plane. 
Therefore, the integral, can be written in the form 

— — J — — — C(s) ds 

2irj (s-b-jw) 


where the path of integration encloses the right half plane. This 
is in the form of Cauchy's integral formula"^ and is evaluated as 
-AC(jto+b) . 

Equation (2-14) thus can be written as 

£a + H(j w)j C ( jw) - A C(jto+b) = G(jw) R(jm) (2-15) 

The problem now is to determine a function C(jw) which satisfies 
Equation (2-15). Unfortunately, a method of determining the 



solution of Equation (2-15) has not been found. The most 
promising approach has been to represent C(joi) as an infinite 
series; however, no real success has been achieved by this 
method. Although the above approach has not yielded an analytical 
solution for the time-varying system, it does come the closest of 
any method attempted. 

An attempt to obtain an analytical expression for the system 
output from either Equation (2-13) or (2-14) has been unsuccessful 
Thus, a different approach must be taken to develop a representa- 
tion which may be used for design purposes. 



CHAPTER III 


DEVELOPMENT OF SYNTHESIS PROCEDURES 
3.1 Introduction 

In the last chapter an attempt to obtain a solution to the 
time-varying system equations in the form of an explicit function 
of frequency proved to be unsuccessful. Design procedures based 
on the knowledge of such functions, therefore, cannot be developed. 
The approach in this chapter is to develop equations and design 
methods which lend themselves to numerical solutions of the 
system equations. This approach should be more successful since 
numerical solutions can be obtained from analog and digital 
computer simulations. The difficulty is the development of 
design procedures which can utilize such solutions. 

A degree of success has been achieved in that two design 
procedures based on numerical solutions have been developed. 
Unfortunately, as was pointed out in Chapter I, a general 
stability criterion is not available for linear time-varying 
systems. Thus, it is not known whether the designs obtained by 
these procedures are stable. When the plant has only a time- 
varying gain, however, sufficient conditions for stability do 
exist so that in such cases it is usually possible to arrive at a 
stable design. 

The two synthesis procedures will be presented in the 
following sections together with illustrating examples. 
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3.2 Synthesis Procedure One 

This design procedure is based on a comparison of the actual 
system step response to a desired system step response. The 
approach is to design the system such that the magnitude of the 
difference in the Fourier transforms of the actual system step 
response and the desired step response falls within some 
specified level. For this method the specifications are assumed 
to be in the form 

I E | = | C-C 0 | < M (3-1) 

where C 0 is the Fourier transform of the desired step response 
and C is the Fourier transform of the actual step response. (See 
Section 1.3.) 

The primary reason for selecting such an approach is that the 
system equation describing | E | lends itself to the development of 
a practical design method. The procedure requires that a series 
of designs be made in which each new design comes closer to 
satisfying the specifications than the previous design. This 
iterative procedure is continued until either a satisfactory 
design is obtained or until the system becomes unstable. 

The system design equations will now be developed. The 
system is shown in Figure 3.1. 
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Figure 3.1 
System under Study 


In order to simplify the notation, let 


X = PY 


(3-2) 


be the symbolical representation between the input Y of a block 
and its corresponding output X. The inverse representation is 
denoted by 

Y = P _1 X. (3-3) 

Using this notation the system equation. Equation ( 2—14 ) » can be 
written 


P _1 C = GR-HC . 


(3-4) 


Let the inverse time-varying plant be written in the form of the 
sum of a time-varying component AP"! and a time- invariant 
component P 0 - ^. That is, 


p-i = p. 


-1 


+ AP -1 . 


(3-5) 
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Also, the system output C can be written as 

C = C Q + E (3-6) 

where C Q is the desired output and E is the error. Substituting 
these expressions for P - ^ and C into Equation (3-4) one can 
ob t ain 

P -1 E + P 0 -1 C Q + AP -1 C 0 + HC 0 + HE = GR. (3-7) 

The prefilter G is now defined such that if the plant becomes P , 
the system step response is the desired response C 0 . Thus, the 
following equation is also satisfied 

P 0 - 1 C Q + HC q = GR. (3-8) 

Employing the relationship expressed in Equation (3-8) , Equation 
(3-7) becomes 

P 0 -1 E + HE = -AP _1 (C 0 +E). (3-9) 

Equation (3-9) is the design equation which is employed in 
this design procedure. The equation is expressed in full 
notational form as 

P 0 _1 (jw) E(j“) + H(jto) E(jto) = 

00 , 

-/ AP - l(jw,jy) [ C 0 (jy) + E(jy) ] dy (3-10) 

— 00 

where P Q ^(jw) is the transfer function of the time -invariant 
portion of the inverse plant. As was discussed in Section (2.4) 
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this equation is a singular integral equation and no general 
method of solution is known. Thus, it is not possible to solve 
for E as an explicit function of frequency. Nonetheless, the 
equation can be used to arrive at a design procedure. 

Equation (3-10) is written in the form 


E ( j to) 


P -1 (jco) + H(jai) -« 


00 r 

/ AP“1 ( joi , jy) C 0 (jy) + E(jy) dy . (3-11) 

—co L 


Applying the specifications given in Inequality (3-1), one obtains 


E (jto) | = 


-1 




/ AP 1 (jto , jy) [" C 0 (jy) + E(jy)l dy 

— 00 -J 


(3-12) 


< M. 


From the above relation it appears that if system stability can 
be maintained, it may be possible to reduce the magnitude of the 
error to any desired level by increasing the magnitude of H. 

Note that H is the only term over which the designer has direct 
control; It is shown in Appendix A that except in one case the 
magnitude of the error can indeed be reduced to as small a value 
as desired by increasing the magnitude of H provided system 
stability can be maintained. The one exception is at a) = o and 
with the plant not having an integrator. That -is, if the plant were 
time-invariant, the system would be classified as type o 

It is possible to design a system which will meet the specifi- 
cations simply by increasing the magnitude of H to some 



arbitrarily large value. This is not a satisfactory approach, 
however, because as the magnitude of H is increased, the trans- 
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mission of the noise from the noise source N to the plant input 
is increased proportionately. Thus, it is desired to increase 
the magnitude of H no more than absolutely necessary in order to 
satisfy the specifications. The design method presented in this 
section will not directly yield the optimal design, but the 
method does provide a systematic approach to increasing the 
magnitude of H. In order to approach the optimal design, the 
designer must utilize the insight gained from the application of 
the design procedure. As has been mentioned, the procedure 
involves an iteration technique where a design is made and then 
tested. If the specifications are not satisfied, a subsequent 
design is made based on the previous design. The subsequent 
design will come closer to satisfying the specification than the 
previous design. 

In order to begin the iteration procedure an initial design 
must first be made. Any stable design could be used as the 
starting design; however, such an approach could yield excessive 
iterations and would not lend itself to the systematic study of 
the system. One possible approach would be to assume that the 
plant were time-invariant but that the parameters could take on 
any constant value within their ranges of variation. The design 
could be made under this assumption using the methods presented 
by Horowitz. 
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Another approach is to begin with Equation (3-13) and assume 
that the integral term involving the error is zero. If the time- 
variations were known as functions of time, the remainder of the 
equation could usually be solved analytically and a design made 

based on this analytical solution. 

The second approach of obtaining a starting design will be 
presented with the design procedure; however, this is not to 
imply that the second approach is superior to the first. The 
second approach is given because it is not included in the 
literature elsewhere as is the first approach. The second 
approach is also used in the examples given in this section 
because the time-variations are assumed to be functions of time. 


The design method will now be presented. For convenience, 
the design equation is rewritten in Equation (3-13) as 


E(jw) 


X(joo) 


CO 

I 0 (ja>) + / AP -1 (ja3,jY)E(jY) dY 
—00 


(3-13) 


where 


X(jw) = - 

P 0 -1 (ja>) + H(jco) 


(3-14) 


and 


00 

I 0 (jw) = / AP _1 (ju,jY)C 0 (jY)dY • 
—00 


(3-15) 



First Iteration: 


The function to be determined is A(jm). For the first 
iteration let A ( joj) be denoted by A 0 (jm). 

Step 1 
Assume that 

CO 

/ AP _1 (ju) J jY)E(jY)dy = 0 (3-16) 

—oo 

so that Equation (3-13) becomes 

E(jw) = A 0 (jto)I o (jw). (3-17) 

Applying the specifications one obtains 

|E(jw)| - |X 0 (jw)| |l 0 (jo))| < M(w) (3-18) 

which by using Equation (3-14) can be rewritten as 


I n (jw) l 1 

— < |P 0 A (jw) + H( jw) | . (3-19) 

M(oj) 


Step 2 

Since the term 1 1 (joj) | depends upon the plant parameter 
variations, it will take on a range of values in the complex 
plane. This range is determined for several frequencies of 
interest. Regions are then determined in the complex plane 
which correspond to the values of I Q ( jco) /M(to) . From these 
regions a region of acceptable values of H(jco) is determined 
in the complex plane for each frequency of interest. 
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Step 3 

Determine an H(jw) as a function of jto such that it lies 
within the acceptable regions found in Step 2 and such that 
| H (j to) j is as small as possible. This step is actually an 
exercise in curve fitting and the H(jm) which is selected 
will depend upon the desired simplicity of the compensation 
network as well as the experience of the designer. Let the 
H(jw) selected in this step be denoted H 0 (jw) . 


Step 4 

The prefilter G(jm) is determined from Equation (3-4). 
Solving Equation (3-4) for G(jrn) one obtains 


G( jm) 



(joi) + H o (jw) 


CpO) 

R(jw) 


(3-20) 


which completes the first design. 

Since H 0 (jm) was designed under the assumption expressed in 
Equation (3-16), it will in all probability not be a satisfactory., 
design. In order to determine whether or not this first design 
satisfies the specifications the design is simulated on an analog 
or digital computer and the frequency domain representation of the 
error is obtained. Let the error of the first design be denoted , 
E 0 (jm). Note that the simulation is merely the means of obtaining 
the solution to the equation 

E Q (jw) = \,(jw) [l Q (jw) + / AP -1 (jw,jy)E Q (jY) dyj . 


(3-21) 
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The magnitude of E Q (jm) can now be compared to the specifi- 
cations M(w) . If E 0 (jw) does not meet the specifications over 
part of the frequency range, another compensation which will be 
denoted AjCjrn) must be designed. In order to see the philosophy 
behind the design of Aj(jai) let the error which corresponds to 
the design be written as 

E 1 < JY) = E 0 (jy) + ejCjy). (3-22) 

The equation which gives EjCjto) can then be written 

Ei(jw) = 

(3-23) 

00 00 — 

+ / AP '’ 1 Cju,jY)E o (jY)d Y + / AP _1 (ju,jY)ei(jY)dY . 

L -» -00 J 

The first integral in Equation (3-23) is obtained from Equation 
(3-21) as 
00 

/ AP _1 (jm,jY)E 0 (jY)dY = - I 0 (Jy). (3-24) 

The value of the second integral in Equation (3-23) is unknown; 
thus, as with the assumption expressed in Equation (3-16) it is 
assumed that 

CO 

/ AP _1 (jm, jY)e 1 (jY)dy = 0. (3-25) 

— 00 ' 

Substituting Equations (3-24) and (3-25) into Equation (3-23) and 
applying the specifications one obtains 
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Ai(joj)l | E Q ( joj) 
1 * 0 (jw) | 


< M(w) . 


(3-26) 


This inequality is then used to determine A^(jm). 

Second Iteration : 

Step 1 

Simulate the design which was obtained in the first iteration 
on an analog or digital computer and obtain the maximum value 
of the magnitude of the frequency domain representation of 

Vv' 

the error which is denoted E om (jto) for the frequencies of 
interest. 

Step 2 

Determine the maximum value of Xj(joi) from the relation 
| A , < jtu) | < | A < jco) | for E om (ju>) > M(co) (3-27a) 

|* 1 (jw)|.< |x 0 (ju)| for E orn (ju)) - M W' (3-27b) 

Employing equation (3-14) this relation can be rewritten as" 


Ufa) S IVhja) + H l«“> I f ° r e ob <» ■> 

° ' 1 (3-28a) 

J A ( < Iv 1 ^^ + H l(j“>l for E om (j«) < M(tt). 

° (3-28b) 
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Step 3 

Determine the region of acceptable H(jw) in the complex plane 
from Inequality (3-28) for the frequencies of interest. 

Step 4 

Determine an H^(jm) as a function of jw such that it lies 
within the regions of acceptable H(jco) found in Step 3 and 
such that |H^(jw)| is as small as possible. 

Step 5 

Determine the prefilter G(jm) from Equation (3-20). 

The subsequent iterations are identical to the second itera- 
tion. If the compensation can be designed to yield a stable 
system at each iteration, the method will eventually yield an 
acceptable design since |X(jm)| of each new design is smaller 
than |x(joj)| of the previous design. The method itself implies 
nothing concerning the stability of the system so that a sta- 
bility criterion must be applied as a separate part of the 
design. Any stability criterion which is expressed in the 
frequency domain can be integrated into the design procedure by 
using the criterion to modify the regions of acceptable H(jco) 
while the basic procedure remains unchanged. 

An example will now be presented. The purpose of the 
example is to illustrate the mechanism of the procedure so the 
design itself will not be carried to completion. Discussion of 
the important considerations in setting specifications and 



49 


i 

in choosing the desired transfer function and the nominal plant 
is also given in this example. 

Design Example One 

Let the plant be described by the first order differential 
equation 

f(t)c(t) = y(t) (3-29) 

where 

f(t) = k x + k 2 [l - e -a(t_T) J u(t-x) (3-30) 

and 

0.1 < k L < 1 

0 . 1 < kj^ + < 1 

0.1 < a < 10. 

It is assumed that the value of t is unknown. From the above 
statement it is seen that the function f(t) can vary exponentially 
between any value in the range from . 1 to 1 with no restrictions 
on when the variation can take place. A block diagram of the 
system is shown in Figure 3.2. 



Figure 3.2 

Block Diagram of System in Example 
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Let the desired step response be given by 
c Q (t) = (l-e -t )u(t) . 

In specifying c D (t) one should insure that it is possible for 
the plant to deliver such a response. That is, the designer 
should insure that there exists a reasonable plant input such 
that the actual output of the plant will be close to the desired 
output. This qualitative check is an obvious one since no 
amount of feedback can make the plant deliver what it is not 
capable of delivering. For this example if the plant output 
c(t) were to equal c Q (t) then the plant input would be given by 

y (t) = f(t) e -t u(t) 

which is a function that imposes no unreasonable demands, in 
theory at least, on the compensations G and H. That is, y(t) 
contains no impulses or derivatives of impulses. 

A second consideration is the specification M. At high 
frequencies the time-varying plant can be approximated by a 
time-invariant plant ^ so that at high frequencies |E(jw)| can be 
rewritten 

<M(u>). (3-31) 

The magnitude of the functions G, R, H, P, and P Q must all go to 
zero as w becomes large and since 1 >> | PH | then 


E = 


PGR 
1 + PH 


_ P 0 GR 

1 + PH 
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e| & |p-P 0 | I GR | < M(w) . 


But by Equation (3-20) GR is chosen as 


GR = (P^ 1 + H)C o 


so that 


JL + H 

?o 


|e| ^ |cp-p 0 ) [ 



(3-32) 


(3-33) 


- 1 + PH - P 0 H < M(u) . (3-34) 

p o 

Note that since P Q is time-invariant then 



Since | PH | << 1 and | P q H | << 1 then one can write 

|E(Jw)|^ - 1 C ( jw) <M(w). (3-35) 

From this formulation it is evident that M(co) , C 0 (jm) and P Q (jw) 
cannot be chosen in a completely arbitrary manner if the specifi- 
cations at high frequency are to be met. Note that if a poor 
selection of M is made such that M goes to zero faster than C Q at 
high frequency, the specifications can never be satisfied at the' 
higher frequencies. This difficulty can be avoided if M is 



chosen as a constant at high frequencies. With M as a constant at 
high frequencies, the specifications can always be met by choosing 
the magnitude of the feedback compensation H large enough to 
satisfy the specifications over the lower frequencies (provided 
stability can be maintained) up to the point where Inequality 
(3-35) will be satisfied by C 0 going to zero. Intuitively it 
appears that there should be no difficulty in letting M be some 
arbitrarily small number at high frequencies although this may 
not be true in all cases. At any rate, Inequality (3-35) must be 
satisfied in the high frequency range. 

The specification for the example must now be chosen. 

Since frequency domain specifications have received little study 
there is no definite basis on which to choose the specifications 
for the example. Rather than attempt to relate time domain and 
frequency domain specifications at this point, specifications 
will be chosen arbitrarily. Let M for the example be chosen as a 
constant of value 0.06. Thus 

M(u>) =0.06 


so that 

|E(ju) | = | C(jio) - C 0 (ju>) | < 0.06 

The nominal plant P 0 (jw) must now be chosen. No definite 
rule can be deduced concerning the selection of P Q (jio); however, 
a guide line can be determined upon examination of Equation (3-13) 
Equation (3-13) can be written in the form 
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E(jto) =. 

CO 

^cjto) / rp-i( JMijY) - p 0 -i (jto,jY) 

— CO L 

The object is to choose |A(jw)| as large as possible and still 
have | E ( jui) | fall within specifications. Obviously, the smaller 
the value of the magnitude of the integral, the larger |x(joi)| 
can be for a given | E (jm) | . The nominal plant P 0 _1 (jw , j y) 
should be chosen such that the magnitude of the integral will be 
as small as possible. 

The design equation for this example will now be determined. 
Note that the differential equation describing the plant is in 
the general form of Equation (2-9) so that the frequency impulse 
response of the inverse plant is evaluated as shown in Equations 
(2-10) and (2-11). The frequency impulse response of the inverse 
plant is found to be 

k _ f 

p-1 (jy,jw) = jw[k 1 + _?.] 6 ( j y - j w ) 



I 




JY - jw jy - jw + a 


The nominal plant must now be determined. Let the transfer 
function of the nominal plant be given by P Q (jw). Referring to - 
Section 2.2, since the nominal plant is time-invariant, the fre- 
quency impulse response of the nominal plant is given by 

Po'^jY.jw) = P Q -1 ( jw) 6 ( jy - Jw) 
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so that AP 1 (jy,jw) is found to be 


AP' 


’kjYjw) = + ^2 J j 


to - P 0 -1 (jto) 


5 C Jy - jto) 


+ jtok 2 e“OY - jto)r 


2tt 


jy - jto jy - jto + a 


(3-36) 


Substituting the above expression for AP -1 (j y,jto) into Equation 
(3-13) and rearranging terms , one obtains 


E(jY) 


[ k i + — ] jy + H(jy)~] = - |"[ 


L 2 J L L 

2 


3Y “ P 0 " 1 (jY) 


c 0 Cjy) 


oo 

+ f k 2 e-(jy - jw)r 
2ir 


1 


jy - jto jy - jto + a 


C ( jto) dio . 


It is desired to select P Q ^(jy) such that the magnitude of the 
right hand side of this equation is as small as possible. The 
only term in the equation that involves P 0 ~^(j<o) is the first 
term on the right hand side. The possible values of the second 
term which is the integral term are not known so that the "best" 
p 0 -1 (j<o) cannot be determined. As an alternative P 0 -1 (jto) will 
be chosen such that the magnitude of the first term will be a 
minimum. Since kj + k 2 /2 can range from .1 to 1, the nominal 
plant is chosen to fall half way between .1 and 1 or 
P o -1 (j“) = 0.55jw. Thus AP -1 (j y, jto) is given by Equation (3-36) 
with P 0 -1 (jio) = 0.55j(o. 

The first step in the first iteration of the design procedure 
is to determine I 0 (jy) which is given by 
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CO 

I Q (jy) = / AP -1 (jY,ju)C 0 (jw)dto 

—CO 


I (jy) = / P -1 (jY, jo))c o (jw)daj - / P 1 (jy,jw)c o (joj)da). (3-37) 

^ —m 


The first integral in Equation (3-37) is the Fourier transform of 
the input to the plant when the plant output is c Q (t) and the 
second integral is the Fourier transform of the input to the 
nominal plant when the output is c Q (t) so that the expression for 
I (jy) can be written as 


I n (3Y) = / 


(kj - 0.55) 


+ k 2 [l - 




e _t u(t) e _ 3Y t dt . 


Evaluating this integral one obtains 


i 0 (jy) = 


M - 0.55 + k f 1 - 1 

jto + 1 Ljw + 1 jio + a + lj 

« 

ki - 0.55 +k r ]_ 1 1 e -(l + ja>)T 

jw + 1 [ju + 1 jm + a + lj 


T<0 


T>0 


(3-38) 


The second step in the first iteration is to determine the 
regions in the complex plane for several frequencies of interest 
such that H (jto) falling within these regions insures that 
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Inequality (3-27) is satisfied. Inequality (3-27) can be 
written as 

|I 0 W| < M(w) = 0.06 

|P 0 -l(ja>) + H(jw) | 

which is equivalent to 

|I°( J “)[ llgWj < Q nfi 

I 1 + P 0 (j“) H (j W ) I 

Letting 


P Q C = L 0 (jw) 


the inequality can be expressed as 


I 0 (ja)l [PpW 

0.06 


< 1 1 + L 0 ( jo>) | . 


The problem now becomes that of determining acceptable regions for 
L 0 (jw). Notice that the acceptable regions of L 0 (jw) take on the 
form of circles in the complex plane centered at the point -1 
and having radius equal to |P 0 (jw)| | I om (jm) | /0.06 where 

is the maximum value of |l 0 (ja>| for the particular frequency of 
interest. Using a digital computer, the maximum value of I 0 (jto) 
was obtained from Equation (3-38) for values of w equal to 0.5, 

1, 2, and 5. The boundaries of the acceptable regions of L Q (jco) 
for these values of L 0 (jw) are shown plotted in Figure 3.3. Note 
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that the acceptable regions of L 0 (jw) lie on the outside of these 
circles. The third step in the design procedure is to determine 
a L 0 (jco) which will just lie within the acceptable region found 
above. The function selected for L Q (jto) is 

L 0 (jm) = 24.8(jm + 2) 

+ 1) (jco + 4) 

and is also shown plotted in Figure 3.3. 

The fourth and final step in the first iteration of the 
design procedure is to determine the prefilter G(jrn) from the 
equation 

G(jm) = [l + L 0 (jm)l CoW . (3-39) 

P 0 (jm)R(jm) 

This is a straightforward matter of substituting L Q (jm) ,P 0 (jto) , 
and C 0 (ja)) /R(jco) into the expression for G(jw) and simplifying. 

The expression for G(j(o) was found to be 

G(jw) = Q-55(ja)) 3 + 2.75(jco) 2 + 18 . 56 (jco) + 32,73 
(jco) 3 + 6 ( j co) 2 + 9(jm) + 4 

In order to determine whether or not this design actually 
meets the specification, the system was simulated on an analog 
computer and | E Q ( j to) | was obtained at selected frequencies for a 
number of step responses using different values of kjyk^, a 
and x . The values of | E Q (jco) | were determined directly from 
the analog computer using the method presented by Dick and 
Wertz^ for obtaining the Fourier transform of a signal on an 
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analog computer. An envelope of the maximum values of | E Q ( jco ) | 
is shown plotted in Figure 3.4. Observe that the specifications 
are not met and that a second iteration must be made. 

The first step in the second iteration which is the simu- 
lation of the design and the obtaining of |E 0 (jm)| has been 
accomplished. The second iteration will start directly with the 
second step. The second step in the procedure is determing the 
values of | + Hj(jw)| such that Inequality (3-28) is 

satisfied for a number of frequencies of interest. Inequality 
(3-28) can be written in the following form 

I 1 + ^o 1 l E p(j U) ) 1 < |l + Li(jco) | for | E q ( jto) | > M(co) 

M(w) 

|l + L 0 (jw)| < |l + L-^ (jw) | for | E q (jrn) | < M(w) 

which, for this example, becomes 

|l+L 0 (jw)| [ E 0 ( j co) | < 1 1 + Li (jw) | for |E ri (jm)| > 0.06 (3-40a) 

0.06 

1+L 0 (ju)| < |l + Li(joi)| for | E q ( jco) | < 0.06. (3-40b) 

Note that the left hand side of the above inequalities can easily 
be determined from the plot of I shown in Figure 3.4 and 

the plot of L q (jco) shown in Figure 3.3. 

The third step is the determination of the regions of 
acceptable values of L^(jm) in the complex plane. These regions 
are found from Inequalities (3-40) and are shown plotted in 
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Figure 3.5 for several values of to. The boundaries take on the 
form of circles centered about the point -1 with radius equal 
to the larger of the two values |l + L Q (jto) | l E om (jw) | /0 .06 or 
|l + L 0 (jto)|. The acceptable regions lie on the outside of these 

circles . 

The fourth step is determining L^jto) such that it lies 

within the acceptable regions found in step three and is as small 

in magnitude as possible. The I^Cjto) which was selected is 

476 (jto + 1.5) [(jto) 2 + 6(jto) + 36] 

Li(jrn) = j(jo(ju) + 1) [(jto) 2 +12(jio) + (12) 2 ](j<o + 12) 

The polar plot of this function is shown in Figure 3.6. 

The fifth and final step in the second iteration of the 
design procedure is determining G(jto) from Equation (3-39) 
which was found to be 

G ( jto) = 

,55(jto) 5 + 13.75(jto) 4 + 433.4(jto) 3 + 3,072(jto) 2 + 12,732jto + 14,138 
(jto) 5 + 26 (jto) 4 + 337 (jto) 3 + 2 ,328 (jto) 2 + 3,744j(o + 1,728 

This completes the second iteration of the design procedure. 

In order to determine whether or not this design meets 
specifications, the system was simulated on an analog computer and 
| E-^ (jto) | was found for a number of step responses corresponding to 
different values of kj, k 2 , a, and t at selected frequencies. The 
envelope of these is shown in Figure 3.7a. It is seen that the 
specifications are satisfied; thus, the design may be considered 
acceptable. 
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Although this completes the first example, some comments 
are in order. Even though the system meets specifications it 
must be noted that the design is not necessarily the best design 
possible. As has been pointed out, the best design is that 
design which will meet specification with as small a value of 
H(joi) as possible in order to minimize the noise transmission 
from the plant output to the plant input. When the best design 
has been obtained the design inequalities will become equalities 
except at high frequencies. In a system where it is particularly 
important to have |H(joi) | as small as possible, more iterations 
may be made using the inequality 

|i + l„(»| |yj»>l < |! +Ll ( 1m) | 

M(co) 

for all values of oj . The subscript zero represents the present 
design and the subscript one represents the subsequent design. 

Using this procedure the designer should be able to gain con- 
siderable insight into a particular system and thus arrive at a 
design which is closer to the best design. It should also be 
noted that system stability was not considered in the example. 

If an applicable stability criterion is determined, the criterion 
will place restrictions on.H(ju) which can then be incorporated 
into the procedure by modifying the acceptable region of H(ju>) or 
equivalently L(jw) in the complex plane. 
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A second design example will now be presented which illus- 
trates a modification that can be made in the design equation 
although the procedure basically remains unchanged. 

Design Example Two 

The plant of this design example is described by the second 
order differential equation 

g(t)x(t) + f(t)x(t) = y (t) 

where 

f(t) - Fj + F 2 (l - e" a(t " T l ) )u(t - t L ) 
g(t) = Gj + G 2 (l - e - ^ ^ “ T 2^)u(t - t 2 ) 

and 


0.01 

< 

G 1 

< i ; 

0.01 

< 

G 1 + G 2 < 1 

0.1 

< 

F 1 

< 0.25 ; 

0.1 

< 

F x + F 2 < 0.25 

0.1 

< 

a 

< 10 ; 

0.1 

< 

b < 10 


It is also assumed that and t 2 are unknown. The desired step 
response is given by 


c Q (t) 


= [jL — e ^ *^ t cos (l ,414t) + 0.98 sin(1.414t)J u(t) 


The frequency domain representation of c o (t) is 


C Q (jw) 


4 

[o u>) 2 + 2.8jw + 4] 
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The specifications are arbitrarily chosen as 
M(o)) = Maxf 


* 

•N 

5 

; 0.01 

(jw + 1) (jw + 5) 





A plot of these specifications for low frequency is shown in 
Figure 3.7b. Note that the specifications are consistent with the 
considerations discussed in Example One. That is, the desired 
plant output does not impose unreasonable demands on the plant 
input and it is possible for | H ( jto) | to go to zero at high 
frequency and still meet the specifications since M(w) goes to a 
constant value at high frequency. Following the procedure 
illustrated in Example One, the nominal plant was chosen as 


P 0 _1 <jw) = 0 . 505( jw) 2 + 0 . 175 jw 


and the frequency impulse response of the inverse plant was 
found to be 


P -1 ( JY ,j<*>) 


(G X + (jw) 2 + (Fj_ + £2 )j u 

2 2 


<5 C JY " jw) 


(jo,) 2 e~tty ~ j“> T 2 + jo)e"<JY - 

2 n (jy - jw) 

(jw) 2 e-(3Y ~ j^T^jp, e -(jy - jw) # (3-4D 

2ir(jy - jw + b) 2tt ( j y - jw + a) 


Let the portion of P -1 (jy,jw) which is multiplied by 6 ( jy - jw) 
be denoted by P a _1 (jw). Note that P a -1 ( jw) 6 ( jy-jw) has the form 
of a time- invariant plant (see Section 2.2), while Pb _1 (jy,jw) 
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represents a time-varying plant. Thus, separating the frequency 
impulse response of the inverse plant into components . 

P ""l(ja))6(jY - jw) and P b - ^ (jy, jw) effectively separates the 
plant into its time-varying and time-invariant components. 

The function AP~^ (jy , jw) can now be written as 

AP -1 (jy,jto) = P -1 (j y ,jm) - P 0 -1 (jy, jw) 
or 


AP -1 (jy,jio) 


p a -1 (j<*>) " p 0 _1 0^) J A ( jy - jto) + P b " 1 (jy,jm). 

(3-42) 



Figure 3.7b 

Low Frequency Specification 
Of Design Example Two 


The design equation, Equation (3-13) , can be written as 
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E(jY) : P -l(j Y ) + H(jy) / AP _ 1 (j Y ,jw) [c Q (jw) + E(joj)] dm . 

— 00 

(3-43) 

Substituting Equation (3-42) into Equation (3-43) one obtains 


-1 


00 -i 

/ J^p-iQo)) - P“ 1 (joi) J 6(jY-jo))E(j(i))dto 


E (jY) = P 0 “i(j Y ) + H(j Y ) 

o 

00 

+ / AP”-*- (ji0 ,j Y) c 0 (j w ) du + / P b - -*-(jY».j<d)E(ju>)du) 


(3-44) 


Evaluating the first itegral and solving for E(jy)> Equation 
(3-44) can be written 


E (j y) = Uiy) 


I 0 (jY) + / P b -1 (jY.jw)E(jw)du 


(3-45) 


where 


I 0 (jy) = / 4P 1 (jY,ju))C 0 (j(jj)do) 


(3-46) 


and 


X(j Y ) = — 

P a -1 (jy) + H(jy) 


(3-47) 


This equation can now be substituted in place of Equation 
(3-13) in the synthesis procedure. The advantage in using 
Equation (3-45) can be seen by comparing the terms which are 
ignored in the first step of the first iteration. By using 
Equation (3-13) the entire integral 
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/ AP 1 (ju,jY)E(jy)dy 

is assumed equal to zero (see Equation (3-16)), while by using 
Equation (3-45) only the term 

00 

/ P b -1 (jy, jm)E(jco)da) 

—00 

is assumed equal to zero. 

The first step in the design procedure is to assume 
00 

/ P b -1 (jY,ju))E(jw)dto = 0 (3-48) 

—CO 

then apply the specifications to E(jy) as follows 
| E( jy) | = | l(jy)I 0 (jy) | < M(y) . 

Using Equation (3-47) this can be written 

I I o(jy) I < lP a -1 (jy) + H(jy) 1 . (3-49) 

M(y) 

Note that both I (jy) and P 0 -1 (jy) depend upon the plant parameters 
so that the ranges in the complex plane of both of these functions 
must be considered in determining the acceptable range of H(jy). 

It must be observed at this point that I c (jy) and P a ’'"(jy) are not 
independent. However, to simplify the task of determining the 
acceptable regions of H(jy) from Inequality (3-49) it will be 
assumed that I Q (jy) and P a ~^(jy) are independent. This assumption 
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should not introduce any more error than is inherently present in 
the first iteration due to the assumption expressed in Equation 
(3-48). 

Since (im) is obtained by inspection from Equation 

(3-41), it only remains to determine I 0 (jy) which is given by 

CO 

I (jy) = / AP -1 (jy» jm)C 0 (jm)dm 


or 


00 

I 0 (jy) = / P" 1 (jY,j^)C 0 (jm)daj 

—CO 


00 

/ P o _1 (jY,jm)C o (j0j)d(i). 

— oo 


0 


This is equivalent to the Fourier transform of the plant input 
with the plant output equal to c Q (t) minus the Fourier transform 
of the nominal plant input with the nominal plant output equal to 
c Q (t). Thus, I 0 (jY) is given by 


= / 


G 3 + G 2 [l 



where 


^4e -1 * 4t cos (1 ,414t) - 3 .96e -1 ,4t sin(l ,414t)J u(t) 
+ [f 3 + F 2 [l - " T l)] u(t - t x )J 

[2.83 e- 1 * 4t sin(1.414t) u(t)]l e'^t ;dt 


G 3 = Gj - 0.505 
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and 

F 3 = ¥ l - 0.175. 

Performing this straightforward but lengthy integration, a compli- 
cated expression for I D (jy) is obtained consisting of four 
functions for different ranges of and ^ 2 ' This expression is 
given in Appendix B . 

The second step in the synthesis procedure is to determine 
the regions in the complex plane at selected frequencies which 
correspond to I D (jy) and P a “^(jm) for different values of the 
plant parameters. From these two regions corresponding to 
I 0 (jY) and P a “^(jy) one can determine acceptable regions of H(jy). 
This was done for y equal to 1, 5, 10, 20, 30, and 50. The 
resulting acceptable regions for H(jy) are shown in Figure 3.8. 

The third step is to select an H(jy) which falls within the 
acceptable regions determined in Step 2. Because of the limited 
size of the analog computer available, it was necessary to limit 
the feedback compensation to third order. The compensation which 
was chosen is 

H(jy) = 31.54 (jy + 0.5) (jy + 5) (jy + 50) 

(jy + 2) (jy + 30) (jy + 80) 

The polar plot of H(jy) is also shown in Figure 3.8. 

The fourth and final step in the first iteration is the 
determination of the prefilter G(jy) from Equation (3-28) which 


was found to be 



100 
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Boundaries of Acceptable Regions of H ( j to) 
with Polar Plot of H(ju;) 
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G(jy) = 

2 (iv) 5 + 227(1 y) 4 + 5497(jy) 3 + 1853(iy) 2 + 38367jy + 15769 
(jy) 5 + 115(jy) 4 + 2938(jy) 3 + 12584(jy) 2 + 23920jy + 19200 

The system was simulated on an analog computer and E(jy) 
was found for a number of step responses corresponding to 
different values of the plant parameters. The envelope of the 
maximum value of these E(jy) is shown plotted in Figures 3.9a 
and 3.9b. Observe that the specifications are met except at 
very low frequencies. This illustrates the basic procedure so a 
second iteration will not be made. Note that |e| is reasonably 
close to the specifications except at low frequency. Thus, the 
second H would be similar to the first except for a higher gain 
at low frequencies. 

Although the designs in the two examples were not carried to 
completion, it appears that there is no difficulty in carrying out 
the steps in design procedure one. It remains to be seen how 
easily a design that is close to optimal can be achieved through 
the procedure., A basic assumption in the procedure is that the 
regions of acceptable L(jw) are circular as shown in Figures 3.3 
and 3.5. In the second example it was possible to arrive at 
acceptable regions which were not circular; however, this was only 
true for the initial design. It would be necessary in subsequent 
designs to assume circular regions. If the regions of acceptable 
L(jw) are, in fact, circular, especially in the region of L(ja)), 
then the design procedure should lead to a design that is close 
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Figure 3.9b 

Envelope of | E(ju>)| for 10 < U) < 20 
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to optinal. If the assumption of circular regions is grossly in 
error, it would be more difficult to arrive at a design which is 
close to optimal. 

The second design procedure will now be presented. Like the 
first procedure, the second procedure is an iterative process. 
However, the assumption that the regions of acceptable L are 
circular is not made in the second procedure so that a design 
which is close to optimal may be able to be obtained in a more 
straightforward manner than is possible in the first procedure. 

3.3 Synthesis Procedure Two 

The procedure is similar to the synthesis procedure presented 
by I. M. Horowitz^ for time- invariant systems having ignorance, 
with the exception that iterations must be performed for a time- 
varying system. In this procedure the specifications are given 
in the form 


: 1 - 


C(jai) 

C 0 (jw) 


- k 2 


(3-50a) 


9'l < arg [c(ju) - C 0 (jw)J < 0 2 . (3-50b) 

The design equation for this procedure will now be developed. 
Referring to Section 2.4, the system is described by the equation 


H(jo))C(ju) + Y(j u ) = G(ja))R(jco) 

where T(j^) is the Fourier transform of the input to the plant. 
Factoring C(jw) from the left hand side of the equation one 


obtains 
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H(j w ) + — ■ 

C(joO 


C(jaj) = G(jcu)R(jw) 


(3-51) 


Under the conditions that the plant is at its nominal time- 
invariant value the system will be described by the equation 


H(jto) + L 


Po (jw) 


C 0 (joi) = G( jm) R(jm) . 


(3-52) 


Dividing Equation (3-51) by Equation (3-52) one obtains 


f H(jco) + ^>1 C(jm) 
J- C(jm)J 

. ,(W + dw] c ° (W 


= i. 


Multiplying the numerator and denominator of the right hand side 
by P 0 (jw) and letting 

P 0 (ju))H(ju)) = L 0 (ju) 

the equation can be written 


C(,jai) 

C 0 (j“) 


1 

Y(,ju) 

C(ju) 


+ L 0 (jq)) 
‘ P 0 (jw) 


+ L 0 (jco) 


(3-53) 


Letting the ratio of the plant output to the plant input be 
denoted by , 


P (jw) = 

6q Y (j «) 
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Equation (3-53) can be written 

C(jm) = 1 + LoW . (3-54) 

C oW SnjO)-- + L (jco) 

P eq (jc) 

For the time-invariant case with ignorance P (jm) is the 
transfer function of the unknown time- invariant plant. Using the 
procedure presented by Horowitz, it is possible to determine 
acceptable regions for L 0 (jm) from the region in the complex plane 
corresponding to the possible values of P eq (jo)). For the time- 
invariant case the values of P e q(jw) are dependent only upon the 
plant parameters so that knowing the ranges of these plant 
parameters it is possible to determine the range of P eq (ju)) in 
the complex plane. However, for the time- varying case P eq (ju)) 
is dependent not only upon the plant and its variations but also 
upon the plant input which in turn depends upon the feedback 
compensation H(jm) . Thus, P e q(jw) depends upon L Q (jw) in the time- 
varying case. This presents difficulty since L 0 (jw) is not known 
prior to the design. 

To circumvent this difficulty an iteration procedure is 
employed where an initial design is made and tested. If the 
specifications are not satisfied, a second design is made based 
on the region of P 0 /P eq corresponding to the first design. If 
additional designs are required, each is based on the regions of 
P o/ p eq corresponding to the previous design. In this scheme the 
designer must use his judgment because there is no guarantee that 
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the procedure will lead to a satisfactory design if mechanically 
followed. Intuitively, however, it appears that the procedure 
should lead to a satisfactory design. If large changes in H are 
required, the procedure will at least call for changes which are 
in the proper direction to improve the design. If relatively 
small changes in H are needed, the regions of P Q /P e q should show 
little change from one design to the next so that regions of P D /P e q 
which are very close to being correct will be used in the design 
with the result that a proper design will be made. It is con- 
ceivable that difficulty could arise in starting with a design 
which requires large modifications in H. There is nothing in the 
procedure which prevents an oscillatory effect. That is, the 
first design may require a large increase in |h| and the second 
design may call for a large decrease in |h| back to the first 
design. Of course, in such a case the designer would apply his 
judgment and not decrease |h| as much as specified. It is 
important that the designer consider all previous designs and 
their corresponding regions of P 0 / p eq w ^ en m aking a new design to 
avoid one that is similar to a design which has previously been 
shown to be unacceptable. 

The steps in the second design procedure will now be 
presented. As in design procedure one, the initial design is a 
starting point for the iteration procedure. As was discussed in 
design procedure one, the initial design may be any stable design; 
however, it would be better to start from a more definite design 
such as one based on a time- invariant plant with ignorance. The 
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approach to obtaining an initial design which is given is based 
on the assumption that the plant output is the desired plant 
output for all variations. It is also assumed that the time- 
variations are known as functions of time so that analytical 
calculations can be made. 

First Iteration 
Step One 
Assuming that 


C(jw) = C 0 (joi) 

determine P 0 (jui)/P e q(jm) from the expression 


P 0 (jai )_ P o (jw)Y(j0)) _ P 0 (ju>) 
p eqOw) C 0 (ju>) C 0 (jw) 


CO 

/ P -1 (jw,jY)C 0 (jY)dY 
— 00 


( 3 - 55 ) 


for a number of frequencies of Interest. Plot the range of 
P Q (jw) /Peq(joj) in the complex plane for the possible plant 
parameter variations. 

Step Two 

Using the procedure outlined by Horowitz^ determine the 
acceptable regions of L Q (ju) at a number of frequencies of 
interest and design a L 0 (jw) which falls within these 
acceptable regions. 

Step Three 

Determine the prefilter G(jw) from the relation 


G(jto) 


] 


1 + L, 



c o(3“) 


P 0 (jw)R(j“) 


( 3 - 56 ) 
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The preliminary design is now complete. 

This first design is simulated on a computer in order to 
determine whether or not the system satisfies its specifications. 

If the preliminary design is unsatisfactory, which presumably it 
will be, a second iteration must be made. 

Second Iteration 
Step One 

From the computer simulation of the previous design, determine 
the range of P Q /P e q in the complex plane for the possible 
plant parameter variations. 

Step Two 

Repeat steps two and three of the First Iteration. 

Subsequent iterations are made until a satisfactory design is 
obtained. It should again be emphasized that the designer must 
consider not only the ranges of P 0 /P e q corresponding to the last 
design, but also the ranges corresponding to all previous designs. 
What the designer must guard against is an iteration leading to a 
design which has previously been found to be unsatisfactory. There 
is nothing inherent in the design procedure itself which prevents 
such an occurrance so the designer must consider all previous 
ranges of ^ Q /^ e q in each new design. An exmmple illustrating this 
synthesis procedure will now be presented. Again, the purpose of 
the design is to illustrate the procedure rather than arrive at an 
optimal design for the particular system considered. 
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Design Example Three 

The plant which was considered in Design Example Two will be 
used in this example. Since a design was made in Design Example 
Two, it will be of interest to compare the regions of 
P o0“)/P e q^ ai ^ with the L 0 (j“) which was obtained in Design 
Example Two. In order to utilize the design of Design Example 
Two, the specifications of that example must be put in a form 
compatible with this design procedure. That is, the specifica- 
tions must be put in the form 


k l * 


c(.m) 


C 0 (jw) 


< M-2 


6 1 - Arg “ C 0 (jm)J< 0, 


(3-5 7a) 


(3-57b) 


instead of the form 


|c(j«) -C 0 (jM>| <M(u>). (3-58) 

Note that Inequality (3-58) defines, a region in the complex plane 
in which C(jto) must fall. Likewise, Inequalities (3-57) also 
define a region in the complex plane in which C(jw) must fall. 

The specifications are given in Figure 3.7 so that one must, 
determine the values of k^, k 2 > 0^, and &2 such that the region 
defined by Inequalities (3-57) corresponds to the region defined 
by Inequality (3-58) . The shapes of the two regions are not the 
same so that values of k^, k£ , 0^, 02 can only be chosen such that 
the region defined by Inequalities (3-57) is approximately the 
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same as the region defined by Inequality (3-58) . The values which 
transfer the specifications of Figure 3.7 to specifications which 


can be 

used with 

Inequalities (3-57) 

are listed in Table 

3.1. 

to 

k l 

k 2 

e i 

e 2 

1 

.36 

1.55 

- 35° 

-35° 

5 

0 

3.6 

-180° 

180° 

10 

0 

11 

-180° 

180° 

20 

0 

24 

-180° 

180° 

100 

0 

2.5 (10) 3 

o 

o 

00 

iH 

1 

180° 

400 

0 

1.6 (10) 7 

-180° 

180° 

500 

0 

3.1 (10) 7 

1 

1—* 

00 

o 

o 

o 

o 

00 

r— H 


Table 3.1 

System Specifications 

The first step of the first iteration is to determine the 
range in the complex plane of the values of P Q ( jto) /P eq (jto) under 
the assumption that C(joi) = C Q (jto). The expression for ; 

P (j(o)/P_,_(jto) is given in Appendix C. The regions in the com- 
plex plane corresponding to the values of P Q ( jto) /P^ (jto) for to 
equal to 1, 10, 20, 100, 400, and 500 are shown in Figures 3.10a 
and 3 . 10b . 

Step two of the procedure is to determine the acceptable 
region for L Q (jto). Using the specifications given in Table 3.1 
the acceptable regions for L Q (jto) were determined from the regions 
shown in Figures 3.10a and 3.10b. These acceptable regions for 
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L q ( joj) are shown in Figures 3.11a and 3.11b. Also shown on the 
figure is the polar plot of L^(jco) which was designed in Design 
Example Two. It is known that this compensation satisfies the 
specifications except at very low frequency; however, it is seen 
from the figure that L^(jto) does not fall within the acceptable 
regions. Two conclusions can be drawn from these figures. The 
first is that at least for this problem the approximation 
C(jw) = C 0 (j w ) is a poor one. The second is that a compensation 
which is designed on the basis of these acceptable regions will in 
all probability more than meet the specifications since its 
magnitude would be considerably larger than the magnitude of 
(jm) • 

A second iteration will now be made. However, instead of 
designing a compensation based on the regions of Figures 3.11a 
and 3.11b the L^(joj) which was designed in Design Example Two 
will be used. The system was simulated on an analog computer and 
C(joj) and Z(joj) were evaluated from this computer simulation for 
a number of variations of g(t) and f(t). The values of 
F 0 (joj)Z(jcu)/C(jco) for these variations were then plotted in the 
complex plane. This plot of points was used to determine the 
range of P Q (jw) / p eq ( jto) . These ranges " of P Q ( jco)/P eq ( jw) were 
determined for values of oj equal to 1, 5, and 10 and are shown in 
Figure 3.12. The range of P 0 (jto) /P e q(jco) for higher frequencies 
could not be determined due to practical limitations of the analog 
computer available. These ranges were next used to determine the 
acceptable regions of L 0 (joj) which are shown in Figure 3.13. 



(mf) L l- / “It 
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Q) 



Figure 3.11a 

Low Frequency Boundaries of Acceptable Regions 
of L(ju») with Comparison of L ( j cjo) 
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Figure 3.13 

Acceptable Regions of L(jcu) 
with Comparison of L^(joj) 
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A polar plot of L-^(jw) is also shown in Figure 3.13. It is 
first noted that L^(jw) evaluated at uj equal to one does not fall 
within the corresponding acceptable region for L 0 (jai). This is 
not surprising since the specifications in the last example were 
not met at low frequency using L^(ju). However, at w equal to 
5 and 10, L^Cjco) does not fall within the corresponding acceptable 
region even though it is known that L^(ju) does meet the specifi- 
cations in Design Example Two at these frequencies. This is 
partially due to the fact that the specifications listed in 
Table 3.1 are not exactly equivalent to the specifications given 
in Figure 3.7; however, the primary error comes from the analog 
computer calculations and in determining the region of 
PqC jw) /P pt1 (jm) from these calculations. 

The next step in the procedure would be to design a second 
L 0 (j(o) based on the regions shown in Figure 3.13. Such a design 
will not be carried out since the basic method has been illustra- 
ted; however, note that the magnitude of L 0 (ju>) of such a second 
design would be greater than the magnitude of L^Cjto). Thus, a 
second design would be a step in the right direction toward 
meeting the specifications as given in Table 3.1. 

3 . 4 Summary 

Two design procedures which rely on the numerical , solution of 
the system equation via computer simulation have been presented. 
Simple design examples have also been carried out which demon- 
strate that it is possible to perform the required steps in the 
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procedures; however, an optimal design has yet to be obtained. 

It is of interest to consider the similarities and differ- 
ences in the two procedures. Even though the procedures appear 
to be quite different they are actually very much alike. Note 
that in both procedures the design of the nominal loop gain L Q 
is carried out in the same manner in that acceptable regions are 
determined and an L 0 is designed to fall within these acceptable 
regions. Also, it is shown in the next chapter that the specifi- 
cations for the two procedures can be made identical. Thus, the 
only basic difference . in the two procedures is that in procedure 
one it is assumed that the shape of the regions of acceptable L Q 
are circular, while in procedure two the shape of the regions of 
acceptable'L Q are determined in the procedure itself. Observe 
that because of this assumption procedure one is more easily 
applied than procedure two. In procedure one it is only necessary 
to determine the maximum value of |e| while in procedure two it is 
necessary to determine the complete region of P Q /P e q» It is seen 
in the next chapter that the price paid for the simpler calcula- 
tion in procedure one is the greater difficulty in arriving at an 
optimal design. 

A major difficulty in the procedures is the determination of 

the maximum value of 1 E I and the region of P /P . The method 

11 o eq 

which is presently employed is to run representations of all 
possible time-variations. Unfortunately, even where there are 
only one or two time-varying parameters a large number of runs 
are required. If a system has several time-varying parameters, 
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the number of combinations of time-variations are prohibitively 
large to test representative samples of all possible time- 
variations. Thus, for systems having several time-varying 
parameters a procedure is needed for determining the worst case 
variations without having to perform a search over all possible 
variations . 

In this chapter the examples demonstrated that it is possible 
to carry out the steps in the procedures using numerical solutions 
from computer simulation. However, the examples do not have the 
depth necessary to observe the advantages and disadvantages of 
the procedures nor do they illustrate how a design close to the 
optimal design may be obtained. In the next chapter the design 
of a system having a time- varying gain is carried out in depth 
so that greater insight into the procedures can be obtained. 



CHAPTER IV 


DESIGN EXAMPLE 

4.1 Introduction 

Previous examples have demonstrated the steps in the two 
design procedures with little consideration given to the practical 
problems encountered in carrying out these steps or of the perform- 
ance of the system involved. In this chapter a design example will 
be considered with greater depth of study in an effort not only to 
gain further insight into the design procedures but also to consider 
the practical problems associated with the procedures, as well as to 
investigate the performance of the system. The example which will 
be considered is shown in Figure 4.1. 



Figure 4.1 
Design Example 

The function f(t) is a time-varying gain which can vary between 1 
and 10. The general structure of f(t) will be discussed in more 
detail later in the chapter. Although this example may appear 
rather restricted, note that it is a special case of a more general 
class of systems where the only time-varying element is a time- 
varying gain as shown in Figure 4.2. 
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Plant 

Input 


Plant 



Figure 4.2 

General Structure of Plant Where Only 
Time-Varying Element Is a Gain 

The portions of the plant Pj and P 2 are linear and time- 
invariant with the time-varying gain represented by f(t). Although 
Pi and P 2 are integral portions of the plant and cannot be modified 
by the designer, any system having such a plant can be put In the 
form shown in Figure 4.3 so that P^ and P 2 can be considered as 
integral portions of the compensations G and L. 



Figure 4.3 

General System with Time-Varying Gain 
In addition since ? 2 (s) is a known transfer function at the 
system output, system specifications on C can be obtained from the 
system specifications on C. The portions of the plant Pj and P 2 , 
therefore, need not be considered in the study of the time-varying 
characteristics of the system. Also, since G(s) is primarily a 
prefilter, the salient properties of the class of systems which 
have only a time-varying gain can be investigated by considering 
only the time-varying gain with a feedback loop for compensation. 
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Results of the design example shown in Figure 4.1 can thus be 
applied to the more general class of systems as shown in Figure 4.3 
where f(t) varies by a factor of ten. The transfer function has 
been included in the plant of the example so that the wave forms of 
the time-invariant system will be familiar and the effects of the 
time-variation will be more clearly evident. 

4.2 Methods of Calculations 

Since the design procedures require numerical calculations of 
the Fourier transform for various signals in the system, it was 
necessary to derive a method for performing these calculations. Due 
to the fact that both design procedures require a search over the 
possible time variations for several frequencies, any method used 
must not only be accurate but also must be able to be performed 
quickly in order that a fairly large number of calculations may be 
made within a reasonable period of time. 

The first approach was to simulate the system on an analog 
computer and calculate the required Fourier transform directly on 
the computer. This technique, which is described in Reference (44^ 
is relatively fast and the results are satisfactory as long as the 
value of the Fourier transform is above the noise level of the 
computer and the run time of the computer can be accurately timed. 

A disadvantage of this approach is that a separate run is required 
for each value of frequency unless a large computer is available 
with at least six amplifiers for each frequency of interest in 
addition to the system simulation requirements. A major difficulty 
in analog computation, however, is the synthesis of time variations. 
The method of producing variations must be of such precision that 
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the variations can be duplicated for repeated calculations of the 
Fourier transform at different frequencies, yet it must be flexible 
enough so that the characteristics of the variations may be quickly 
altered. The ability to precisely control and quickly change 
exponential variations on the analog computer is the primary reason 
for these type variations in the examples of the previous chapter. 

The computer system used for this first approach consisted of a 

Beckman 2200 analog computer and an SDS model 920 digital computer 
which had been combined into a hybrid system. The Beckman analog 
computer is a fully shielded, solid state, ± 100 volt machine. The 
particular configuration consisted of a total of 36 amplifiers and 
four multipliers along with a complement of logic and an accurate 
real time clock. In addition, there were sufficient servo potentio- 
meters to provide satisfactory hybrid operation. The SDS digital 
computer had an 8K storage with a high speed paper tape reader and 
teletypewriter input-output. 

Use of the analog computer alone proved to be, on the whole, 
unsatisfactory. It is required that the operator make each run, and 
read and record the values of the calculated Fourier transform. 

This in itself is a lengthy procedure; but, in addition, calcula- 
tions must be made on each data point and the result plotted. The 
entire process proved to be extremely lengthy. 

By utilizing the full hybrid capabilities of the computer, 
calculation times were significantly reduced. The digital computer 
was programmed to change the time variations, control the analog 
computer, read the computed values, perform the required calculation, 
and print the results. Thus, many more calculations could be 
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performed with the hybrid configuration in a given time than with 
the analog alone. 

The analog computer's method of calculating the numerical 
Fourier transform proved satisfactory for the first design procedure; 
however, in the second design procedure difficulties were encoun- 
tered with this method of computation. In the second design pro- 
cedure it is necessary to compute the Fourier transform of the 
system output. Since the system input is a step function, the 
system output generally goes to a constant value which will be 
denoted C M . In such a case the numerical computation methods yield 
the sum of the Fourier transform and the constant C ro . Thus, in 
order to determine the Fourier transform it is necessary to subtract 
the constant Q*, from the computed result. For the majority of 
frequencies it was found that the difference in and the computed 
results was small. This means that the inaccuracies, which are 
small in comparison to the large value of the computed result, are 
large when compared to the small value of the difference in and 

the computed result. Thus, in many cases it was impossible to 
obtain an accurate value of the Fourier transform of the system 
output with the analog method of computation. 

In an effort to eliminate this problem with the second design 
procedure, it was decided to use the digital computer to calculate 
the numerical Fourier transform instead of the analog computer. 

This approach would allow for a more accurate computation of the 
transform. The procedure was to simulate the system on the analog 
computer, sample and store the sample values of the signals of 
interest, which in the case of. the second design procedure are the 
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plant input and plant output, then compute the values of the Fourier 
transform at the frequencies of interest. The algorithm described 
in Reference (48) was used to compute the numerical Fourier trans- 
form of the sampled signals. At this time the system was simulated 
on an EAI-TR20 analog computer with a PDP9-L digital computer being 
used to sample the signals and perform the numerical calculations. 
The combination could not be considered a hybrid computer since the 
digital computer could not be used to control the analog computer. 
Although a number of calculations were performed on this set of 
computers, the attempt on the whole was considered a failure. The 
lack of success was not due to the methods of computation which were 
basically sound but to the inadequacy of the computers themselves. 
This attempt illustrates the fact that relatively sophisticated 
computers are required for the necessary computations in the two 
design procedures . 

The necessary simulations and Fourier transform calculations of 
the example in this chapter were done wholly on a digital computer. 
The computer used was an IBM 1130 with a core storage capacity of 
16K with three disk drives. The basic program used to simulate the 
system was the Continuous System Modeling Program for the IBM 1130 
computer. The program was modified by the addition of several 
subroutines for the purposes of storing signals, calculating the 
Fourier transform, modifying the time-variations , and printing 
results. The primary disadvantage to this approach is that it 
requires a good deal more time to make a calculation than with the 
hybrid computer. To illustrate the difference in run time, assume 
it is necessary to calculate the Fourier transform at five fre- 
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quencies for a signal of ten seconds duration. On a hybrid computer 
with the Fourier transform calculated on the analog portion of the 
computer, one run would require ten seconds of computation time with 
five runs being required. Add approximately 30 seconds additional 
time to allow for data reading and parameter changing and the total 
time becomes 80 seconds. The digital computer requires from two to 
four minutes to generate the required signals of a third or fourth 
order system plus approximately 30 seconds to calculate the Fourier 
transform at each frequency for a total run time of approximately 
five minutes. This time increases substantially as the order of the 
system increases. The primary advantage in using the digital 
simulation is the increase in accuracy that is obtainable. However, 
one also gains other advantages. Time-variations are easily 
simulated and modified on the digital computer so that it is 
possible to develop search techniques for determining the worst case 
variations which is required in the first design procedure and for 
generating randomly varying variations which is required in the 
second design procedure. Thus, once the program is started the 
operator need do nothing more until a solution is obtained. 

4.3 Specification 

As was discussed in Chapter One, system specifications are 
more commonly given in the time domain than in the frequency domain; 
thus, there is no reservoir of experience upon which to draw in 
assigning frequency domain specifications as there is for time 
domain specifications. In order to determine realistic frequency 
domain specifications for the example considered in this chapter, it 
will be necessary to first consider the acceptable system step 
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response in the time domain and then attempt to relate these to the 
frequency domain. The approach will not be rigorous but will 
consist of examining the magnitude of the Fourier transform of a 
sampling of error functions and assigning a frequency domain 
specification based upon this examination. 

Consider the first form of frequency domain specification 

which is given by 


| C(joj) - C (jto)| = | E( jrn) | < M(m) 


(4-1) 


where C Q is the desired step response, C is the actual step 
response, and M is a real number which can depend upon u • The form 
of this specification naturally suggests a time domain specification 

of the form 

|c(t) - c 0 (t)| = |e(t)| < m(t). 


There are two graphical interpretations of this relation. The 
specification given by 

I c(t) - c 0 (t) | < m(t) 


(4-2) 


corresponds to a symmetrical envelope about the desired step 
response within which the actual step response must fall. Such a 
specification is shown in Figure 4.4. 
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Time Domain Specifications as Envelope 
About Desired Step Response 


The specification can also be expressed in the form 

|e(t)| < m(t) (4-3) 

which corresponds to a symmetrical envelope about the time-axis as 
shown in Figure 4.5. 



Figure 4.5 

Time Domain Specifications as Upper 
Bound on Error 


In this form the specification represents an upper bound on the 
error. Although the specifications given in inequalities (4-2) and 
(4-3) are equivalent, the specification as given in inequality (4-3) 
is somewhat more convenient to study since it is independent of the 
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desired s tep response. In all probability, however, the specifica- 
tion for any particular system will first be established in the form 
of inequality (4-2) since one is generally more familiar with 
acceptable step responses than with acceptable error functions. 

For the example considered in this chapter, it will be assumed 
that the desired system transfer function is given by 

T(s) = 

This , transfer function will occur in the system of the example if 
the gain is constant at 4 and the compensations are given by 

H(s) ’ ' Htt 

The transfer function corresponds to a second order system having a 
damping ratio of 0.7 and an undamped natural frequency of 2. The 
step response of this transfer function, which is the desired step 
response, is shown in Figure 4.6. 

In assigning the specification envelope, consideration must be 
given to the usual time domain specifications of rise time, over- 
shoot and settling time. The rise time determines how quickly the 
lower bound of the envelope must increase while the peak overshoot 
determines the maximum values of the envelope's upper bound and the 
settling time determines how quickly the envelope must converge to 
the final value. In the final analysis, however, it will be neces- 
sary to assign an acceptable minimum step response and an acceptable 
maximum step response. These maximum and minimum step rfesponses will 

constitute the specification envelope. Note that the time domain 
specification given in inequality (4-2) demands that the envelope be 
symmetrical about the desired step response. If the maximum and 
minimum step responses are not symmetrical about the desired step 
response, the desired step response must be redefined such that this 
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symmetry is achieved. The assumed specification envelope for the 
example is shown in Figure 4.6 with the desired step response. 

This envelope corresponds to a minimum rise time of approximately 
1.5 seconds, a peak overshoot of 30% and a maximum settling time 
of 8.5 seconds. 

The upper bound on the error which corresponds to the specifi- 
cation in Figure 4.6 is shown in Figure 4.7. The problem now is to 
determine the frequency domain equivalent of this specification in 
the form of an upper bound on the magnitude of the Fourier trans- 
form of the error. This frequency domain specification will be in 
a form as shown in Figure 4.8. 



Form of Frequency Domain Specification 
The desired objective is to determine the frequency domain envelope 
M(o>) such that, if the magnitude of the Fourier transform of the 
error is less than or equal to the maximum value M(co) , the error 
will fall within the time domain envelope given in Figure 4.7. 

A difficulty is now encountered which again emphasizes the 
complexity of the problem of obtaining frequency domain' specifica- 
tions from time domain specifications and vice-versa. Assume that 
an error e(t) satisfies the time domain specification of Figure 4.7 
and the magnitude of the Fourier transform |E(ju))| satisfies the 
frequency domain specification of Figure 4.8. Let the error e(t) 
be shifted in time by the amount t such that e(t - x) now lies 
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outside the time domain specification. The Fourier transform of 
e(t-t) is E(jto)e~j UT which has the same magnitude as the Fourier 
transform of e(t) and thus will still satisfy the frequency domain 
specification. The reason for this difficulty is due to the fact 
that e(t-t) is nonminimum phase so that the uniqueness between 
phase and magnitude does not hold. Thus there is no unique rela- 
tionship between the error e(t) and the magnitude of its Fourier 
transform; It is easily seen then that if the frequency domain 
specification is given in the form shown in Figure 4.8 and the time 
domain specification is given in the form shown in Figure 4.7, there 
will always be functions which will satisfy the frequency domain 
specifications and yet violate the time domain specifications. 

In stable systems it is generally the case that the maximum 
error will occur relatively close to the origin, then will settle 
to zero. Otherwise, the actual step response would closely follow 
the desired step response through its maximum change; then, after 
attaining the final value would show a deviation. It would be 
necessary for a system to possess some form of large time delay to 
exhibit such behavior. Since there is a large class of systems 
which do not possess such time delays and therefore will not 
exhibit a delayed error, it will still be assumed that the time 
domain specification is given in the form of Figure 4.7. 

A number of possible error functions and their Fourier trans- 
form magnitudes will now be examined in order to deter mi ne the 
error function characteristics in the frequency domain. 

Figure 4.9a shows four functions which, for different values 
of time, lie on but do not exceed the error bound m(t). Figure 4.9b 
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Figure 4.9b 

Fourier Transform Magnitudes of 
Error Function Series One 
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shows a plot of the magnitude of their Fourier transforms. Also 
shown for comparison is the magnitude of the Fourier transform of 
the time domain envelope m(t). An error function that is exactly 
equal to m(t) has the largest zero frequency component allowable 
since this component is equal to the area under the curve. But 
since m(t) has few variations, its higher frequency component would 
be smaller than might be allowed. Functions 2 and 3 shown in 
Figure 4.9a cross the zero axis leaving about as much negative area 
as positive area and thus would be expected to have a relatively 
small low frequency component. Note that the type of variation for 
all four functions is similar so they would be expected to have 
similar higher frequency components. This is shown in Figure 4.9b 
where all four functions exhibit peaking around uj = 1.5, then 
decrease rather quickly toward zero. Function 2 exhibits the 
greatest amount of variations in the time domain which shows up in 
the frequency domain as a little larger high frequency component 
than is observed in the other three functions. 

Figure 4.10a shows a series of error functions which has a 
relatively large variation then settles fairly quickly to zero. 

The frequency domain representations are shown in Figure 4.10b. As 
expected, all exhibit a relatively small zero frequency component 
since the areas under the curves are relatively small. The higher 
frequency components depend upon the particular variations. If the 
variations show an oscillation effect as does Function 5, the higher 
frequency components will be larger; whereas, if the error deviates 
from zero only once then returns to zero as does Function 6, the 
higher frequency component will be smaller. 
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Shown in Figure 4„lla is a series of functions having relative- 
ly low frequency oscillatory tendencies but remaining within the 
time domain error bounds. The magnitudes of their Fourier trans- 
forms are shown in Figure 4.11b. The oscillatory tendency is 
evident in large frequency components in the region of w equal from 
2 to 3. Note that Function 10 is smaller than the other three, and 
thus its Fourier transform is also smaller in magnitude. 

Figure 4.12a shows two functions which fall within the time 
domain error bound but are highly oscillatory. The magnitudes of 
their Fourier transforms are shown in Figure 4.12b. As expected, 
the oscillations show up as large high frequency components. 

It would also be of interest to observe the effects in the 
frequency domain of a single fast variation. Figure 4.13a shows 
such a function with its frequency domain representation shown in 
Figure 4.13b. 

Although the error functions shown in Figures 4.9 through 4=13 
represent a small sample of the possible error functions, they were 
selected from the functions investigated because they give a good 
representation of the error functions possessing Fourier transforms 
with large magnitudes which are the types of functions that will 
establish the error bound in the frequency domain. The low fre- 
quency bound, that is, in the region from 0 to 1 radian per second, 
will be established by the Fourier transform of m(t) which is shown 
in Figure 4.9b. This is because as was previously mentioned m (t) 
has the largest low frequency component allowable. To establish the 
high frequency bound, it is observed that the magnitude of the 
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Figure 4.11a 

Error Function Series Three 
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Figure 4,1 lb 

Fourier Transform Magnitudes 
of Error Function Series Three 
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Figure 4.12a 

Error Function Series Four 
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Figure 4.12b 

Fourier Transform Magnitudes 
of Error Function Series Four 
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Figure 4.13b 

Fourier Transform Magnitude of 
Function with Single Fast Variation 
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Figure 4.14 

Possible Frequency Domain Specification Envelopes 
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functions shown is less than .03 for u greater than seven radians 
per second unless the functions are highly oscillatory as shown in 
Figure 4.10a. Such oscillatory functions are generally undesirable 
even though they do not exceed the time domain error bound, and 
thus the high frequency bound should be set small enough to exclude 
these undesirable functions. Nevertheless, a high frequency 
oscillation of small magnitude is often considered acceptable; thus, 
instead of assigning a very small value for high frequencies, a 
relatively large value will be selected for the high frequency 
bound. The minimum value of the frequency domain bound at high 
frequencies will be selected as 0.05. 

The frequency domain specification must now be established in 
the mid-frequency range. A minimum frequency domain envelope would 
be established by the magnitude of the Fourier transform of the time 
domain envelope m(t) which is shown in Figure 4.14. Also shown in 
Figure 4.14 are additional envelopes, the inclusion of which will 
allow additional types of functions. The envelope labeled A will 
allow Function 7 of Figures 4.10 and the Function shown in Figures 
4.13. Note the similarity of these two functions. The envelope 
labeled B will allow the functions shown in Figures 4.9 and the 
functions shown in Figures 4.10 with the exception of Function 5. 
Observe that the large oscillatory nature of Function 5 excludes it 
from this set. The envelope labeled C allows all functions shown in 
Figures 4.10 and 4.11. The envelope labeled D is included in 
Figure 4.14 to demonstrate the large magnitude required in the 
region of ui equal to five to six to include the highly oscillatory 
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functions shown in Figures 4.12. The error functions allowed by 
the envelope labeled B will be assumed to be acceptable in the 
example; thus, the magnitude labeled B will be selected as the 
frequency domain specifications. This selection will exclude the 
oscillatory type of functions shown in Figures 4.11 and 4.12 but 
will include functions having a high degree of variations. 

As was discussed earlier, it is possible for the frequency 
domain specifications to be satisfied and the time domain specifi- 
cations to be violated. With this in mind, it would be of interest 
to examine the magnitudes of the Fourier transform of a sampling of 
functions which violate the time domain specification but may be 
expected to satisfy the frequency domain specifications. Such 
functions could not have large oscillations for they would surely 
violate the frequency domain specification at high frequency, nor 
could such functions have large deviation without being oscillatory 
for they would then violate the low frequency specifications. One 
possible type of function that may violate the time domain specifi- 
cation but satisfy the frequency domain specification would increase 
slowly to a value outside the time domain specification then return 
to zero before the frequency domain specifications are violated at 
low frequency. Three possible functions are shown in Figure 4.15a 
with their Fourier transform magnitudes shown in Figure 4.15b. Even 
though functions 15 and 16 violate the time domain specification 
only a relatively small amount, the frequency domain specifications 
are nonetheless violated at low frequency. Function 17, on the 
other hand, does not violate the frequency domain specifications 
while the time domain specifications are violated. This, of course, 
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Figure 4.15b 

Fourier Transform Magnitudes of Functions Which Violate Time Domain Specifications; Series 
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is because the magnitude of function 17 does not become large 
enough to have too large an area under the curve but does not 

return to zero quickly enough to avoid violating the time domain 

specifications. Functions such as function 17 cannot be eliminated 
by altering the frequency domain specifications without also 
excluding the desirable functions. However, the slowly settling 
component in such functions is probably due to poles near the 
origin as would occur if L(jm) (see Figure 4.3) were assigned a 

zero near the origin and the magnitude of L(joj) were not sufficient- 

ly large in the low frequency range. If functions of this type 
occur, the cause will probably be evident and can thus be eliminated 
from the design. 

Another type of function which could possibly violate the time 
domain specification but not the frequency domain specification 
would quickly reach a point outside the time domain specification 
then return to zero to avoid violating the frequency domain specifi- 
cation at low frequency. Two such functions are shown in Figure 
4 16a with their Fourier transform magnitudes shown in Figure 4 16b. 
Even though function 18 violates the time domain specifications, it 
would in fact satisfy these specifications if it were shifted 
properly in time. Function 18 violates the frequency domain 
specification only slightly at u) equal to five radians per second 
while function 19 completely satisfies the frequency domain 
specification. Note that the magnitude of the Fourier transform of 
both functions remains relatively large at frequencies around six 
to nine radians per second. This is undoubtedly due to the fast 
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Figure 4.16b 

Fourier Transform Magnitudes of Functions 
Which Violate Time Domain Specifications; Series II 



rate of variations of the functions . It is likely that such 
functions can be controlled by reducing the high frequency specifi- 
cation. A reduction at high frequency will not be made on the 
frequency domain specification which has been chosen, but careful 
attention will be paid as to whether or not such error functions 
become present in' the design of the example. 

To summarize, the frequency domain specification which is 
shown in Figure 4. 17 has been chosen as a counterpart to the time 
domain specification shown in Figure 4.7. The selection was based 
on a small but what is felt to be a representative sampling of error 
functions. It was also pointed out that an error function satisfy- 
ing the frequency domain specification will not necessarily satisfy 
time domain specifications as shown in Figure 4.15; however, such 
instances would have to be due to some form of time delay within the 
system or a slowly settling component within the system. In such 
instances the particular portion of the system which is causing the 
difficulty would have to be examined separately to see if it can be 
eliminated or its effect reduced to an acceptable level. 

The frequency domain specifications used in the second design 
procedure must now be determined. These specifications have been 
assumed to be given in the form 


ki(oi) < 


C(.jai) 

C 0 (ju) 


< k2(io) 
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Figure 4.17 

Frequency Domain Specification 
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In order to compare the two design procedures, it will be desirable 
to have the specifications as given in Inequality (4-4) to be 
equivalent to the specifications shown in Figure 4.17. Such an 
equivalence is easily established if the two sets of specifications 
are properly interpreted. Note that the specifications given in 
Inequality (4-4) define a region in the complex plane which is 
illustrated in Figure 4.18. 



Figure 4.18 

Region Defined by Second Form of 
Frequency Domain Specifications 


Also observe that the specification given in Figure 4.17 defines a 
region in the complex plane. This region is shown in Figure 4.19. 
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Figure 4.19 

Region Defined by First Form of 
Frequency Domain Specifications 

It is evident that since the regions shown in Figure 4.18 and 4.19 
are basically different, one region can only be an approximation to 
the other. However, the procedure for using the specification given 
in Inequality (4-4) is a graphical procedure and the region shown in 
Figure 4.19 can be used as satisfactorily as the region shown in 
Figure 4.18. Thus, the specifications used in the second design 
procedure will be the regions defined by the frequency domain 
specification given in Figure 4.17 and illustrated in Figure 4.19 in 
place of the region shown in Figure 4.18. 

4.4 First Design 

In design procedures number one and two, the first design is a 
starting point from which a final design is reached. In the 
previous chapter the methods of obtaining a first design assumed 
that the time variations were known as an explicit function in order 
that the necessary integrals could be evaluated. In the example 
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given in this chapter the time variations are not assumed known as 
an explicit function, so that some other method of determining the 
first design must be used rather than that previously given. It 
would be possible to assume any stable compensation as the starting 
design; however, not only may this lead to an excessive number of 
designs in the procedure, but such an indiscriminate approach does 
not lend itself to an orderly study of time-varying systems. The 
most logical starting design would be one based on a time- invariant 
system with plant ignorance. Since the plant is time-invariant, 
such a design can be made analytically; and, in addition, the 
effects of the time variations on the system can be directly 
observed while the additional feedback required to compensate for 
these time variations is clearly evident. 

The first step in the time invariant design is to find the 
acceptable regions in the complex plane for H(jto). However, 
instead of determining the acceptable regions of H(jco) what will 
actually be found are the acceptable regions for 
-L 0 (jw) = -H(jw)Po(jw) where P 0 (ju>) is the nominal plant. In 
designing the time invariant system one should bear in mind that 
the design equations for both design procedure one and two stem 
from the same basic equation. The difference in the two procedures 
is that procedure one bases the design upon the difference in the 
actual system output and the desired system output, while procedure 
two bases the design upon the ratio of the actual system output to 
the desired system output. Since the specifications found in the 
previous section for the two procedures are exactly equivalent, 
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the acceptable regions for -L 0 (jw) found by either procedure one 
or two will be identical. For illustration purposes the acceptable 
regions will be found using both procedure one and two. 

The acceptable regions for -L 0 (jm) will first be obtained for 
design procedure one. The design equation for design procedure 
one is given as 


E(jw) 


-1 

P 0 (j“>) + H(joj) 


AP 1 (jw > jy)C(jY)dy. 

. 


Since the system is assumed time invariant, the design equation can 
be written 


E (joj) = -rP-l(jui) - P'^jm) ] [ Q,(ja)> + E(jo)>] . 

Po _1 (j“) + H(ju) 


Solving for E(jw) and rearranging one can obtain 

1 ^ 1 


E(J .) . y-w c d (». 

p(j U ) + 


(4-5) 


In selecting the nominal plant, observe that | E ( j co [ may be reduced 
by either increasing the magnitude of the denominator of Equation 
(4-5) or by decreasing the magnitude of its numerator. It is the 
design objective, of course, to increase the magnitude of the 
denominator no more than necessary while maintaining as small a 
value of [ H ( j to ) [ as possible. Thus, P D (jto) should be delected such 
that the magnitude of the numerator will be reduced to as small a 
value as possible. The inverse plant is given by 
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1_ _ joi(^a) V + 2) 

P(ju) = k 

where 1 < k < 10 

and the inverse nominal plant is given by 


1 = juli^ + 2) . 

P 0 (jm) k 0 


The expression for the magnitude of the numerator of Equation (4-5) 
becomes 


1 1 


1 to ( .1 to + 2) 

i to ( i to + 2) 



ko 

k 



1 1 



n 

ko ” k 


jw ( jto + 2) 


For any m the maximum value of this expression will be minimized 
for all k when l/ ko is given by 



1. J. 

kmin + kmax 
2 


= ill = °- 55 
2 


or ko = 1.818, 

Thus, the nominal plant is selected as 

1.818 . 

P 0 (j^) = jw(ju) + 2) 

The regions of acceptable -L 0 (jw) can now be found. Rearrang- 
ing Equation (4-5) as shown on the following page 
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1 - 


E(j“) - p (ju>) 


Pp(ja) 

P(jM) 


c oCj“) 


p (jcu) - 1 + 1 + H(joj)p o (joj) 


one can write 


E(ji») = 


- Q( j tn) 


Q(ju) + L 0 (ju) 


C n (j“) 


where 


P 0 Ow) 

Q(j u ) - p (joj) “ 1 


(4-6) 


and 

L 0 (j w ) = H(jO))P 0 (ju) + 1 = L 0 (jw) + 1. 

Applying the specifications to Equation (4-6) , one obtains the 
inequality 

| Q(jw) | 

I Co(J»> I [ Q(» + £ 0 (ju» I i M(») 


or 

| Q(j0)) + L o (j0)) | > I Q(j Ui ) 1 1 C o (j0i) I (4-7) 

M(w) 

At any specific frequency <*> the function Q(jw) takes ori a range of 
possible values in the complex plane depending upon the plant P(ju)). 
For each value of Q(joj) the right hand term of Inequality (4-7) 
takes on a specific value which then specifies a forbidden region 
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for -L 0 (jco). Determining the forbidden regions for all values of 
Q(jco) then determines the forbidden region of -L 0 (jco) at frequency 
co. The compliment of the forbidden region automatically gives the 
acceptable region of -L 0 (jco). To illustrate the procedure, the 
acceptable region for -L Q (jco) at co equal to two will be found. 

The range of Q(jco) must first be determined. The expression for 
Q(jco) is given by 

Q(j(o)= _ x 

P (jw) 

= 1.818 ,i co ( .1 co + 2) _ 1 

jco(jco + 2) k 

= 1.818 

k ~ A * 

Since the constant k takes on all values between one and ten, 

Q(jco) takes on all real values between -0.818 and 0.818. From 
Figure 4.17 the value of M(2) is found to be 0.4 and the value of 
|C 0 (j2)| is calculated as 0.352. Table 4-1 gives the values of 
the right side of Inequality (4-7) for different values of Q(jco) 
at co equal to 2 . 

1 Q( jco) |. 1 C n ( jco) | 

Q M(co) 


-0.818 

0.731 

-0.5 

0.446 

-0.2 

0.179 

0 

0 

0.2 

0.179 

0.5 

0.446 

0.818 

0.731 


Table 4-1 

Values of the Right Side of Inequality (4-7) 
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To determine the forbidden region of -L 0 (jw) for a particular Q, 
consider the value of Q equal to -CL 81 8. At this value of Q one 
finds from Table 4-1 that the magnitude of Q(jw) + L 0 (ju>) must be 
greater than or equal to 0.731. The region where this is true is 
illustrated graphically in Figure 4.20. The distance from Q to 
-L 0 is the magnitude of Q + L 0 , so that the forbidden region of 
-L 0 for Q equal to -0.818 corresponds to the region inside a 
circle with center at -0.818 and radius 0.731. 



Figure 4.20 

Illustration Showing the Determination of 
Acceptable Region of -L 0 (jw) 

Repeating this procedure for several values of Q, one is able to 
determine the forbidden regions corresponding to each value which, 
when superimposed upon one another, gives the total forbidden 
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region. The acceptable region is then the region outside the 
forbidden region. Figure 4.21 shows the total acceptable region 
of -L 0 for (D equal to 2. 



Figure 4.21^ 

Acceptable Region of -L 0 for w = 2 v 

Having the acceptable regions of -L 0 , it is a simple matter to 
determine the acceptable region of -L 0 . Since L 0 is given by 

L 0 = 1 + Lq , 

then 

-L o = 1 “ ^o* 

Thus by adding one to each point in the acceptable region of -L 0 , 
the acceptable region of -L 0 is obtained. This corresponds to a 
shifting of the acceptable regions of -L 0 to the right by one to 
obtain the acceptable regions of -L 0 . 
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The frequencies of interest for which the acceptable regions 
are to be found for the example must now be selected. In calcu- 
lating the Fourier transform on the computer, the frequencies 
which may be calculated are multiples of the fundamental frequency 
given by 

% = ii 

T 

where T is the total system time of the calculation . It will be 
assumed that the system will be close to steady state in ten 
seconds so that T is taken as ten seconds. If, in practice the 
system does not reach steady state in ten seconds, T can be 
doubled or quadrupled until steady state is reached and the 
appropriate harmonic selected. With T equal to ten seconds, to 0 
is equal to 0.628. The multiples of this frequency which were 
chosen are 0, 1, 2, 3, 7, 10 and 15. This selection gives fre- 
quencies fairly equally spaced on a Bode plot and should also 
give adequate representation of the functions of interest. The 
frequencies corresponding to these multiples are 0, 0.628, 1.25, 
1.88, 4.39, 6.28, and 9.42. The acceptable regions of -L 0 (ju>) 
which correspond to these frequencies are shown in Figure 4.22. 

An acceptable region for w = 0 cannot be found since L 0 (ju>) is 
indeterminate at to = 0. If it is found that the error is too 
large at u = 0, a correction will be made in the second design. 

The reduction of the magnitude of L 0 (jto) is not as critical at low 
frequencies as it is at high frequencies since at low frequencies 
the loop gain is adequate to limit the noise transmission to the 
plant input. 
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Figure 4.22 
Acceptable Regions of 
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The acceptable regions for -L 0 (joj) will now be found for 
design procedure two. The design equation for design procedure 
two is given by 


C(i ai) 
C o (j«0 


1 + Lr.CiOl) 


Pq(JM ) 

P(ju) 


+ L 0 (j“) 


The nominal plant is chosen the same as was selected for design 
procedure one. 

P 0 (j“) = 1 -‘ 818 

jm(jo) + 2) 


Since the system is assumed to be time-invariant, the equivalent 
plant is given by 


p eq <j“> 


k 

jw( ja) + 2) 


where 


1 < k < 10. 


The ratio of the nominal plant to the equivalent plant is thus 
found to be 

1.818 

PqW = ju(iu + 2) = 1.818 

p eq(j“) ]E k 

ju(jw + 2) 


so that the design equation becomes 


C 


(jw) 1 + L 0 (ju) 


(j“> 


1.818 


+ 


(4-8) 
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The objective is to chose an L 0 (jid) such that the value of 

C(jw)/C 0 (j«) as given by Equation (4-8) falls within the 

( 28 ) 

appropriate specification regions for all values of k. 

As discussed in the previous section, one form of the 
specification region is a circle in the complex plane of radius 
M( w ) and centered. at the point C 0 (jw). The specifications will 
be satisfied if (C(j^) falls within this region. Dividing every 
point in the region by C 0 (jw), one obtains a region into which 
C(jw)/C 0 (jw) must fall for the specification to be satisfied. 

This region consists of a circle centered at the point one with 
radius M(u>)/ |c 0 ( j ui) |. As with design procedure one, the frequen- 
cies of interest are selected as 0, 0.628, 1.25, 1.88, 4.39, 6.28, 
and 9.42. The regions of acceptable C( j^)/C Q (j o>) for these 
frequencies are shown in Figure 4.23. 

The regions shown in Figure 4.23 are still not in the most 
convenient form for determining the acceptable region of -L 0 (j w ), 
however. To determine the appropriate form, consider Figure 4.24 
as shown on page 144. 

Shown in Figure 4.24 is a hypothetical region of P 0 /Peq and 
a line a . passing through the point one. To determine whether or 
not any point A on the line a lies within the acceptable region of 
-L 0 , one draws the vector from point A to one which is equal to 
1 + L 0 . A vector from A to any point within the regions of P 0 /P e q 
gives P 0 /P e q + L 0 . These two vectors give the numerator and 
denominator of Equation (4-8) so that C/C 0 can be determined. If 
C/C Q as determined by Equation (4-8) falls within the acceptable 
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Figure 4.24 

Determining a Point on the Boundary of the 
Acceptable Region for -L 0 (jw) 

region of C/C Q for all the points within the region of P 0 /Peq> 
point A then falls within the acceptable region of -L 0 . Point A 
lies on the boundary of the acceptable region of -L 0 if C/C 0 just 
falls within the acceptable region of C/C 0 for all P 0 /P e q» The 
object is to determine the boundary points for a number of lines 
passing through the point one. These boundary points then 
determine the region of acceptable -L 0 . 

The determination of the boundary points by the procedure 
just described is a trial and error method and can become fairly 
lengthy. However, if the specification region for C/C D were 
modified so that the numerator appears to be a constant and the 
denominator a variable, which corresponds to the form of Equation 
(4-8) , the determination of the boundary of acceptable -L 0 becomes 



a simpler matter. To illustrate, Figure 4.25 shows the region 


° f p o/ p eq> 


and superimposed is an assumed modified region of C/C, 



Figure 4.25 

Illustration of the Use of the Modified 
Region of C/C Q 

The point one is common to both regions and the vector from 
point A to one which is 1 + L 0 also corresponds to C. Observe 
that if point A is moved away from one along line a , the region 
of C/C Q becomes proportionally larger while retaining the same 
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4. 27a 

:(jw)/C 0 (jo>) : to = 0.628 




150 



Figure 4.27c 

Modified Region of C(ju) /c o (jto): w = 1.88 
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Figure 4.27d 

Modified Region of C(jto) /C D (jto) : to = 4.39 
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Figure 4.27f 

Modified Region of C( ju) /C Q ( jw) : u> = 9.42 
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Figure 4.28 

Determination of a Point on the Boundary 
of the Region of Acceptable -L 0 (jo)) : to = 1.25 
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of acceptable -L 0 ( j u ) . The acceptable regions of -L 0 (j w ) 
obtained by this method are identical to the regions obtained for 
design procedure one which are shown in Figure 4.22. Figures 4.27a 
through 4.27f show the modified regions which correspond to the 
specification regions given in Figure 4.23. 

Since the acceptable regions of -L 0 (jaj) are known, a suitable 
L 0 ( j “) can now be designed for the time-invariant system. It is 
desired to have -L 0 ( j m) just lie on the boundaries of acceptable 
-L 0 at the point where the distance to the origin is a minimum. 

For the acceptable regions shown in Figure 4.22, the points closest 
to the origin lie on the negative real axis. However, for the 
magnitude of LgCjm) to decrease there must be some degree of lag 
associated with L 0 (j<u). Thus, the design of L 0 will be a 
compromise between maintaining as small a magnitude as possible 
while obtaining the necessary rate of decrease. The design of 
was made by trial and error. A preliminary design was 
first made using the Bode plot, then the final design was reached 
by adjusting the preliminary design so that the polar plot of 
-L 0 (J“) came as close as possible to lying on the boundaries of 
the acceptable regions of -L 0 (jw). The expression for the L 0 (s) 
which was selected is 

L ( s ) = 24.17 (s + 0.9) (s + 6) 

s(s + 1.4) (s + 7) (s + 12) 

and has a polar plot as shown in Figure 4.29. The expression of 
H(s) is easily found from the relation 
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H(s) 


LqCs) 

V®) 


and is found to be 


H(s) = 13.3 s 3 + 8.9s 2 + 19.2s +10.8 

s 3 + 20.4s 2 + 110.6s + 117.6 


Generally, the number of poles over zeroes will be greater than in 
the above L 0 (s) and H(s). The high frequency poles were not 
included in order to keep the solution time for the computer 

simulations as short as possible. The high frequency poles will 

be added after a suitable time- varying design has been made. 

The prefilter G(s) is determined from the expression 

G(s) = 1 + L o( s > T 0 (s) 

Po(s) 

Making the appropriate substitutions one obtains 

G(s) _ 2.2 s 5 + 22.4s 4 + 175.6s 3 4- 553.9s 2 + 699,3s + 261.1 

s 5 + 23.2s 4 + 171.7s 3 + 508.9s 2 + 771.7s + 470.4 

Note that the far off poles of L 0 (s) will also be in G(s); however, 
the far off poles in G(s) are not necessary as they are in L(s). 

In fact, the far off poles of G(s) can be assumed to be in T Q (s). 
This will give a simplier expression for G(s) and will result in 
little change in the system response. 

To determine whether or not the design is adequate for the 
system under time-varying conditions, a simulation was made on an 
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IBM 1130 digital computer using the Continuous System Modeling 
Program for the 1130 in conjunction with the algorithm for 
calculating the numerical Fourier transform as given in Reference 
(48). A more detailed discussion of the simulation procedure is 

given in Appendix D. 

The results of the search for the maximum value of the error 
will first be presented. A series of runs was made with the 
function representing the time-varying gain having the general 
form as shown in Figure 4.30. The function has 11 equally spaced 
points starting from time equal zero to nine seconds. The points 
can take on any value between one and ten with the function 
varying linearly between points. Thus, the maximum slope for the 
function is ten except at the selected points. 



Figure 4. 30 

Type of Function for Time-Varying Gain 


In searching for the function which would give the maximum error 
it was quickly found that the larger errors occur when the 
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specified points take on either maximum or mini m um values. The 
variations having the greatest effect on the error occur near time 
equal zero, and any variation after five or six seconds has little 
effect on the error. By far the largest error occurs when the 
time variations start at ten and decrease quickly to one. In 
fact, variations starting at values other than near ten produce 
errors no larger than those observed when the gain is held 
constant at ten, no matter what the subsequent variation. 

It was found that a significant increase in error occurred 
only at frequencies zero and 0.628, while the design appeared to 
be adequate for all higher frequencies. Figures 4.32 and 4.33 
show the step responses and error functions corresponding to 
these maximum errors. Figure 4.37 shows the magnitude of the 
Fourier transform of these error functions. Note that the error 
functions in both cases increase quickly, then slowly decrease to 
the final value with the result of an excess area under the curve. 
What apparently happens in the system during these particular 
time-variations is that when the gain starts at ten, the system 
has the ability to accelerate quickly to a relatively high 
velocity; however, after the system has attained the velocity the 
gain is suddenly reduced to one, so that the rapid acceleration 
effect is lost, thus allowing the system to coast to a higher 
value than would be the case if the gain remained at ten. With 
the gain at the lower value, the system returns to the final 
value more slowly than if the gain were at ten. The effect is 
also observed in Figure 4.33 when the system decreases to zero 
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Figure 4.31 

Functions Resulting in Maximum Error 
During First Series 



Step Response 
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Figure 4.32 

Step Response and Error Function Corres 
To Maximum Error at u = 0 During First 


Step Response 


162 



Figure 4.33 

Step Response and Error Function Corresponding 
To Maximum -Error at u = 0.628 During First Series 
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Figure 4.34 

Functions Resulting in Maximum Error 



Step Response 
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Figure 4.35 

Step Response and Error Function Corresponding 
To Maximum Error at u = 0 During Second Series 


Step Response 



Figure 4.36 

Step Response and Error Function Corresponding 
To Maximum Error at to = 0 During Second Series 
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more quickly after four seconds as a result of an increase in 
gain between four and five seconds. 

In an effort to observe the effect of an increase in slope on 
the error, a second series of time-varying runs was made with the 
same time-varying gain as shown in Figure 4.30, except the 
variations are completed in 4.5 seconds instead of nine seconds, 
thus increasing the maximum gain from ten to 20. A slight increase 
in the error at the two frequencies zero and 0.628 was observed; 
but, again, the design proved to be adequate for the higher 
frequencies. Figure 4.34 shows the variations causing the maximum 
error at the two frequencies and Figures 4.35 and 4.36 show the • 
corresponding step responses and error functions. The magnitude 
of the Fourier transforms of the error functions is shown in 
Figure 4.37. Note the close similarity of these functions to 
those shown in Figure 4.31 through 4.33. 

A final effort was made to cause errors at the higher fre- 
quencies by having the time-varying gain make discrete jumps 
between the values at each of the 11 points instead of varying 
linearly between points. Only negligibly higher values of error 
were observed, however, over those obtained in the second series 
of runs. 

It is felt that an excellent representation of possible time- 
variations were obtained in these runs. A number of preliminary 
runs were first made and analyzed to determine the types of 
variations causing the larger errors. It was these preliminary- 
runs that showed the larger errors to occur when the specified 
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Figure 4.37 

Maximum Error Due to Time Variations 
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points of the function take on either the maximum or minimum 
values and that the largest errors occur when the function started 
at ten and quickly decreased to one. Special functions such as 
oscillations of various frequencies and step functions were then 
selected and run. Finally, a large number of functions were 
selected at random to help insure that some possible variation 
resulting in larger errors had not been overlooked. None of the 
errors due to these special or random selections, however, caused 
errors larger than those observed in the preliminary runs. In 
all, more than five hundred functions were run requiring approxi- 
mately 70 hours of computer time on the IBM 1130. 

The determination of the regions of P 0 /P 6 q for m = 6.28 will 
now be discussed. In reality, the determination of the regions 
of Po/Peq were carried out simultaneously with the search for the 
maximum | E | since both P Q /P e q and [ E | were calculated for 
each function. The procedure for determining the region of 
P 0 /P e q was to plot the values of P G /P e q for each time-varying 

function from which the general outline of the region could be 

) 

observed. It must be recognized that the procedure is not exact 
but, nonetheless, yields useable results. Because of the rela- 
tively large number and wide diversity of functions which were 
run, it is felt that a reasonably good representation of the 
regions of P Q /Peq was obtained. The regions are shown in Figures 
4.38 through 4.43, along with the specification regions of 
acceptable P Q /P e q for the L 0 of the design. The only frequency 
at which the specifications are not met is 0.628, which is in 
agreement with the results obtained in the maximum |e| 
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Figure 4.38 

Region of P 0 /P ea for co = 0.628 





171 



Figure 4.40 

Region of P 0 /P ea for w = 1.88 
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Figure 4.42 

Region of P 0 /P 6 q for a) = 6.28 
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Figure 4.43 

Region of P c /P e q for oj = 9.42 
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determination . 

The regions of P Q /P e q for frequencies 0.628, 1.25, 1.88, and 
9.43 are well concentrated about the region of P 0 /P eq corresponding 
to the time-invariant case. This concentration is similar to the 
observations made in the previous chapter and is no surprise; 
however, the regions of P 0 /Peq for frequencies 4.39 and 6.28 do 
not show this concentration which was an unexpected result. 

The region of P 0 /P eq for 6.28 is especially interesting since 
it extends into the left half of the complex plane and a good 
distance into the lower half plane. In fact, points were found 
as low as -j6. However, such large magnitudes of P 0 /P e q are not 
unacceptable as is seen in Figure 4.42. Note what large values of 
P 0 /P eq imply about the system output C(ju) . Since P Q is a constant, 
P 0 /P eq becomes large only if P eq becomes small. Since P eq is 
defined as the plant output divided by the plant input, P eq becomes 
small only if the time-varying plant does not effectively transfer 
the signal at the frequency of interest. In other words, P Q /P eq 
will be large if the frequency content of the output is small 
while that of the input is not so small. In the example the 
frequency content of the output at 6.28 cannot only be small or 
even zero, but can actually be at some finite magnitude 180° out 
of phase with the desired plant output, thus allowing portions of 
P 0 /P e q to He an the left half plane. 

Figures 4.38 through 4.43 clearly illustrate why the design 
is satisfactory for the time-varying gain at frequencies of 1.25 
and above. The time-variations do have a significant effect on 
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the system as shown by the increased size of the region of P 0 /Peq - 
However, when the specification region of acceptable P 0 / p eq is 
made large enough to enclose the region of P 0 /Peq for the time- 
invariant system, it also encloses a great deal of additional 
area which is large enough to allow for the increased size of the 
region of P 0 /Peq in th e time-varying system. That is, the design 
is not only acceptable for the range of P Q /P e q in the time 
invariant system, but it is also large enough to accept the much 
larger region of P 0 / p e q in the time-varying system. It is 
speculated that if the region of P Q /P e q for the time-invariant 
system came closer to filling the specification region of Po/Pgq 
such as would occur if a plant pole could take on a range of 
values, the region of P Q /Peq for the time-varying system would 
not be able to lie entirely within the acceptable region of P 0 / p e q* 
It is also of interest to observe some of the step responses 
and error functions of the system for time- variations other than 
those resulting in the maximum value of |e|. Seven typical step 
responses with their corresponding time-variations and error 
functions are shown in Figures 4.44 through 4.50. The time 
duration of the randomly selected error functions was placed at 
seven seconds since, as has been previously observed, faster 
time variations and variations after seven seconds have little 
effect on the error. The step responses are well behaved despite 
the time-varying gain. The maximum overshoot appears to be not 
over 30% and the settling time of the function is not over seven 
or eight seconds in the worst case. Although fast rise time is 
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Step Response of Time-Varying System 
Time-Invariant Design 
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Figure 4,47 

Step Response of Time-Varying System 
Time-Invari a nt Design 
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Step Response of Time-varying System 
Time-invariant Design 
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Figure 4.49 

Step Response of Time-varying System 
Time-Invariant Design 
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Figure 4„ 50 

Step Response of Time-varying System 
Time-invariant Design 
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usually not considered an undesirable property, it is fast rise 
time that is responsible for the large values of the error func- 
tions shown in Figures 4.44, 4.46, and 4.49. The time-varying 
gain in the responses having large values of error start at or 
near ten, x^hile functions starting at or near one such as shown in 
Figures 4.45, 4.48, and 4.50 result in step responses having 
little or no overshoot with a relatively slow rise time and thus 
smaller values of error. It is evident from these figures that 
the form of the step response and thus the error function is 
highly dependent upon the starting value of the time- varying gain 
and is less sensitive to subsequent variations. This basically 
is the same observation made earlier in the determination of the 
maximum value of |e|. 

Since the time-invariant system design is unsatisfactory for 
the system under time-varying conditions, a second design must be 
made. The second design will first be made using design procedure 
one . 

4.5 Second Design: Procedure One 

Equations (3-28) are the design equations for the subsequent 
designs required in procedure one. The design equations can be 

for E om > M (4-9a) 

for E om < M (4-9b) 

where M is the frequency domain specification, E om is the maximum 


written as 


J om 

M 


1 + L„ 


1 + L r 


1 + L 


1 + L 
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error of the previous design, L Q is the compensation of the 
previous design and L' Q is the compensation to be designed. 

The minimum values of |l + L^ Q | are now found from Inequali- 
ties (4-9). At zero frequency the magnitude of L' Q is infinite 
so that Inequalities (4-9) cannot be applied. However, dividing 
Inequality (4-9a) by P 0 and setting to equal to zero one obtains 

|H(o)| < |H'(o)| 

which sets a minimum value on |h'(o)|. Table 4-2 gives a 
summary of the specification M, the values of | 1 + L 0 | and the 
corresponding E om of the previous design, and the minimum values 


of |1 

+ L ' Q | for 

the frequencies 

of interest. 


to 

M. 

^om 

U + L ol 

Min | 1 + L ' Q | 

0 

1.13 

1.32 

| H ( o ) | = 1.22 

| H '( o ) | > 1.43 

0.628 

0.65 

0.808 

2.34 

2.91 

1.25 

0.51 

0.543 

1.613 

1.72 

1.88 

0.42 

0.42 

1.34 

1.34 

4.39 

0.11 

0.11 

1.01 

1.01 

6.28 

0.05 

0.05 

0.95 

0.95 

9.41 

0.05 

0.23 

0.93 

0.93 


Table 4-2 

Summary of Specifications 
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The minimum values of |l + L" Q | correspond to regions in the 
complex plane of acceptable -L' Q . These regions are found to lie 
outside circles about the point one with radius equal to the 
minimum value of jl + 1,' Q \ . 

The regions of acceptable -L ' Q corresponding to the values 
given in Table 4-2 are shown in Figure 4.51 together with the 
regions of acceptable -L 0 for the time-invariant system. A 
shortcoming of design procedure one is now evident. Since the 
region of unacceptable -L ' 0 is assumed to be a circle centered 
about the point one, it is possible to exclude -L" 0 from large 
areas in the complex plane which would actually be acceptable, 
thus resulting in an overdesigned system. This will be especially 
true at high frequency where a small increase in the radius of the 
circle results in a much larger required increase in |l' 0 |. 
Unfortunately, it is at these higher frequencies that | L " 0 1 must 
be decreased to as small a value as possible in order to limit 
the noise transmission to the plant input. For example, consider 
the regions corresponding to w equal 6.28. If the boundary of 
acceptable -L 0 for the time -invariant system were not known, the 
boundary of acceptable -L ' Q allows |l q | to become as small as can 
be designed. Such a design could not satisfy even the time- 
invariant system, much less the time-varying system, so that the 
radius of the region would have to be increased which would 
demand a much larger increase in | L "* 0 1 than would really be 
required. 

It was decided to design -L ' Q so that it would lie within the 
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acceptable region of -L^ 0 at to equal 0.628 but would lie approxi- 
mately half way between the boundaries for -L 0 and -L 0 at to 
equal to 1.25 and 1.88. At to equal to 4.39, 6.28, and 9.42 only 
the acceptable regions of -L Q would be considered. The reason for 
this is that the specifications were basically satisfied with L Q 
except at to = 0.628 while a design based solely on the acceptable 
region of -L' 0 would result in an underdesign at the higher 
frequencies. The design which was chosen has an L^ 0 (s) given by 


L ' 0 (s) 


26.88 


(s + 2.2) (s + 5.8) 
s(s +2) (s + 7) (s + 14) 


The polar plot of -L' a is shown in Figure 4.51. The corresponding 
H(s) is given by 

H(s) = 14.78 s 3 + 10s 2 + 28.76s + 25.52 
s 3 + 23s 2 + 140s + 196 

and G(s) is found to be 

G( s ) =22 s 5 + 25s 4 + 212.88s 3 + 744.8s 2 + 1165.07s + 685.98 
s 5 + 25.8s 4 + 208.4s 3 + 680s 2 + 1108.8s + 784 

The higher frequency poles in L/ 0 (s) and H(s) will be added in 
the final design. The above H(s) satisfies the specification on 
H'(o) given in Table 4-2. More lag is required in L' Q than in 
L 0 because of the higher rate of decrease required between the 
frequencies 0.628 and 1.25. 

A comment should be made concerning stability. Since the 
example has only a time-varying gain, sufficient conditions for 
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stability are known. A convenient criterion is the circle cri- 
terion which is given in Reference (38). The system of this 
example will be stable if the polar plot of -l,' does not lie 
within or enclose a circle of radius 0.8181 and centered at point 
one. The above design satisfies this criterion so that the 
system is assured to be stable. 

The new design was tested by simulating the system and 
determining the maximum magnitude of error. The duration of the 
time-varying gain as illustrated in Figure 4.30 was set at seven 
seconds. Since L ' Q is similar to L Q , faster variations should 
have little effect on the magnitude of the error as well as 
variations after seven seconds. As with the first design, it was 
found that the larger errors occurred for variations that start at 
ten and decrease quickly to one and that variations starting at 
values other than near ten produced errors no larger than those 
observed with the gain held constant at ten. It was also verified 
that variations occurring after five or six seconds have little 
effect on the error. 

A plot of the maximum error corresponding to this design is 
shown in Figure 4.52 along with the specification envelope and the 
m^xi.nuim error corresponding to the first design. An important 
observation is that the increase in |l + L'J at a particular 
frequency appears to have substantially affected the maximum 
error only in the vicinity of that frequency. For example, an 
increase of |L Q | at 0.628 did not cause an increase in [ E [ at 
1.88 as could have occurred in a time-varying system since the 
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Figure 4.52 
Plot of Maximum [ E [ 

For Second Design Using Procedure One 


frequencies are not independent. As was suspected, the system is 
overdesigned; that is, |e| is considerably smaller than the 
specifications in the low frequency region. 

The time variations resulting in the maximum error at 10 = 0 
and m = 0.628 are shown in Figures 4.53 and 4.54, respectively, 
with the corresponding step response and error function. Note 
that these variations are almost identical to the time variations 
found for the first design shown in Figures 4.31 and 4.34. The 
step responses and error functions in the second design are also 
similar to those in the first design. In the second design, 
however, the peak overshoot is reduced to 19% from 32% in the 
first design. The settling time in the second design has also 
been significantly reduced. 

Although the design satisfies the frequency domain specifica- 
tion, the error function does not fall within the time domain 
specification envelope. Nevertheless, the step response is well 
behaved and would undoubtedly be considered acceptable. In 
reality, the time domain specification envelope should probably 
have been chosen to allow larger errors in the range 0 to 1 
second since fast rise time is generally not considered an 
undesirable property in a system. To give an indication of the 
step response for other variations, variations are shown in 
Figures 4.55 and 4.56 with their corresponding step responses 
and error functions. 




Figure 4.53 

Time Variations Resulting in Maximum | E j 
0 With Corresponding Step Response and Error Function 




Figure 4.54 

Time Variations Resulting in Maximum [ E [ 

0.628 with Corresponding Step Response and Error Function 




Figure 4.55 

Sample Time Variation with 
Corresponding Step Response and Error Function 
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4.6 Second Design: Procedure Two 

In design procedure two the approach is simply to use the 
regions of P 0 /P eq from the first design to determine regions of 
acceptable -L ' Q for the second design. The regions of P 0 /P e q 
for the first design are shown in Figures 4.38 through 4.43, 
and the specification regions are given in Figures 4.27a through 
4.27f. Using techniques previously described in the first design, 
the acceptable regions of -L 0 are determined as shown in Figure 
4.57. Except at the frequency 0.628, the acceptable regions of 
_L ' Q do not require the magnitude of -L 0 to be substantially 
greater than was required in the first design. This is expected 
since the first design satisfied the specifications at frequencies 
1 .25 and higher. 

A design satisfying the specification regions was made by trial 
and error and the resulting L G (s) f- s gf ven ky 


L V s ) 


,25.55 


s + 5 

s ( s + 7) (s + 12) 


The polar plot of L ' 0 is shown in Figure 4.57. The design 
satisfies the circle criterion for stability which is given in 
Reference (38) so that system stability is assured. The corre- 
sponding expression for H(s) is 

H (s) = 15.15 s 2 + 7s + 10 
s 2 + 19s + 84 


and the expression for G(s) is 

s 4 + 21s 3 + 149.55s 2 + 360.86s + 275.52 


G(s) = 2.2 


s 4 + 21.8s 3 + 141.2s 2 + 311.2s + 336 
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Figure 4.57 

Boundaries of Acceptable Regions of -L' with 
Polar Plot of — Lq ( j to) 
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The high frequency poles will be added in the final design. This 
design requires more lag than the first design because of the 
higher rate of decrease between the frequencies 0.628 and 1.25. 

It should be pointed out that the magnitude of E is somewhat 
smaller in this design than in the second design using procedure 
one. The design was tested by simulating the system and determ- 
ing the regions of P Q /P eq . The duration of the time-varying gain 
illustrated in Figure 4.30 was again set at seven seconds. 
Simulations were run for over two hundred variations with the 
variations being both specifically selected and chosen at random. 
The resulting regions of P 0 /P eq are shown in Figure 4.58 through 
4.63, along with the acceptable regions of P 0 /P eq corresponding 
to L ' Q . The regions of P Q /P eq fall within the acceptable regions 
of P 0 /P eq which means that the design satisfies the specifications. 
However, since the regions of P D /P eq do not touch the boundaries 
of acceptable P D /P eq at all frequencies, the system is overdesigned 
although the overdesign is not appreciable. 

Observe that the regions of P 0 /P e q f° r this second design are 
almost identical to the regions of P D /P eq found for the first 
design. This result is extremely promising for design procedure 
two since, for this example at least, the regions of P D /P eq do not 
change appreciably for changes in L 0 . Thus, designing an L 0 
based upon regions of P 0 /P eq corresponding to a previous design 
does not appear to be an unreasonable approach. 

To give an indication of the type of step responses which 
this time-varying system exhibits, a number of step responses with 
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Figure 4.62 

Region of P 0 /P e a f° r 00 = 6.28 
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Figure 4.63 

Region of P 0 /P e q for ta = 9.42 
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the corresponding error functions and time-variations are shown 
in Figures 4.64 through 4.71. For comparison, the first two 
responses correspond to time— invariant conditions with the gain 
set at ten and one, respectively. As with the earlier designs, 
the time-variations cause an increase in overshoot and settling 
time. The maximum overshoot in the time- varying case is approxi- 
mately 24%, which is an increase over the 19% observed for the 
second design using procedure one and a decrease over the 32% 
observed for the first design. Such a result is expected since 
the magnitude of L' Q for this design is smaller than the magnitude 
of L' Q in the second design using procedure one and is larger 
than the magnitude of L' Q in the first design. 

4.7 Third Design 

In making the second design it was found that a design based 
on procedure one, which is the technique using the maximum |e|, 
will generally result in an overdesigned system because the 
acceptable regions of -L 0 are assumed to be circular. Use of the 
second procedure, which is based on the regions of P 0 /P e q> results 
in the more accurate shape of the acceptable region of -L 0 5 
however, the size of the acceptable region cannot be as precisely 
determined as can be done by the method in procedure one. The 
third design will be based on a combination of procedures one and 
two. The shape of the acceptable regions of -L 0 will be assumed 
to be as was found in the second design using procedure two, while 
the size will be determined from the maximum error. 

The maximum error magnitude was found for the second design 




Figure 4.64 

Step Response of System with Gain 
Constant at 10, Second Design, Procedure 
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Figure 4.65 

Step Response of System with Gain 
Constant at 1, Second Design, Procedure 
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Figure 4.71 

Sample Step Response of Time-varying System 
Second Design, Procedure Two 
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based on procedure two and is shown in Figure 4.72. Also shown in 
Figure 4.72 are the maximum error magnitudes for the second design 
based on procedure one and the maximum error magnitudes for the 
first design. The first design is not satisfactory, while both 
of the second designs result in an overdesigned system. However, 
in the second design based on procedure two the magnitude of L Q 
needs to be decreased only a small amount at the frequencies 0.628, 
1.23 and 1.88 for the error to just lie on the specification 
curve. To determine the amount of decrease needed at to = 0.628, 
it will be assumed that the magnitude of the error varies linearly 
with the magnitude of L Q between the first design and the second 
design based on procedure two. The difference in the magnitude 
of the error of the second design can be increased by .025 which 
is a 16% change based on the 0.156 difference. The acceptable 
region of -L 0 at w = 0.628 for the third design is then determined 
as follows. 

Figure 4.73 shows a plot of the acceptable region of -L 0 at 
0.628 for the first design and for the second design based on 
procedure two in the area of interest. The magnitude of L 0 for 
the first design at 0.628 is 2.02 and for the second design is 
2.61. The difference in the two magnitudes is 0.59. The 
boundary of -L 0 at oj equal 0.628 for the third design will be 
placed at 16% of the difference between these two magnitudes as 
shown in Figure 4.73. 

At w equal 1.25 the magnitude of acceptable -L Q needs to be 
changed only a little; thus, -L Q for the third design will be 
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Figure 4.72 

Comparison of Maximum |e| for 
First and Second Designs 
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chosen to fall approximately half way between the magnitudes of -L Q 
for the first design and -L 0 for the second design based on 
procedure two. At the frequencies 1.88 and higher the magnitude 
of -L will be chosen the same as the magnitude of -L 0 of the 
first design since in the first design the specifications are 
just met at these frequencies. The resulting acceptable regions 
of -L 0 for the third design are shown in Figure 4.73. 

The polar plot of the third design is shown in Figure 4.73 
together .with the acceptable regions of -L 0 . The expression for 
L (s) of the third design is 


L 0 (s) 


24.61 


s + 5 

s(s +8) (s + 10) 


This design satisfies the circle stability criterion so that 
system stability is assured. The corresponding expression for 
H(s) is 

H(s) = 13.54 s 2 + 7s + 10 
s 2 + 18s + 80 


and for G(s) i s 


N n s 4 + 20s 3 + 140.6s 2 + 332.3s + 246.1 
G(s) = 2.2 — • 

s 4 + 20.8s 3 + 134.4s 2 + 296s + 320 


The high frequency poles will be added in the final design. Note 
that this design is very similar to the second design based on 
procedure two. The major difference is the small reduction in the 
magnitude of L 0 over all frequencies. 

The design was tested by simulating the system and determining 
the regions of P 0 /P e q as well as the magnitude of the error. As 
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with the second design, the duration of the time-varying gain 
which is illustrated in Figure 4.30 was set at seven seconds. 

Almost 300 variations were run with variations being both specifi- 
cally selected and randomly selected. It is again felt that 
reasonably accurate results were obtained. 

The maximum error magnitude is shown in Figure 4.74 and the 
regions of P 0 /Pgq are shown in Figure 4.75 through 4.80. As was 
observed in the first design, the maximum errors occurred for time- 
variations starting at ten and quickly decreasing to one. The 
only significant increase in error due to the time-variations 
occurred at the frequencies zero and 0.628. The time variations 
causing the maximum error magnitude at w equal to zero and 0.628 
are shown in Figures 4.81 and 4.82, respectively, along with the 
corresponding step response and error function. Note in Figure 
4.74 that the magnitude of the error is slightly larger than the 
specifications in the region of 0.628. Otherwise, the system 
just satisfies the specifications. Figure 4.75 also shows that 
the specifications are not exactly met at 0.628 since the region 


of 0 ) does not fully lie within the acceptable region of P q /P 


eq* 


Since this system is very close to meeting the specifications 


exactly, it is very likely that this design would be considered 
acceptable. A series of time-variations with the corresponding 
step responses and errors are shown in Figures 4.83 through 4.88. 
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Figure 4.74 

Maximum [ E ] for Third Design 
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Figure 4.75 

Region of P 0 /P eq for “ = 0.628, Third Design 
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Figure 4.77 


Region of P Q /P e q for co = 1.88, Third Design 
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Figure 4.78 

Region of P 0 / p e q ^ or w “ 4.39, Third Design 








Step Response 
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Figure 4.81 

Time-variation Causing Maximum |e| at 




Step Response 



Figure 4.82 

Time-variation Causing Maximum |e 
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Figure 4.84 

Sample Step Response of Time-varying System 
Third Design 





Figure 4.85 

Sample Step Response of Time-varying System 
Third Design 
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Sample Step Response of Time-varying System 
Third Design 
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Figure 4.88 

Sample Step Response of Time-varying System 
Third Design 
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Note that the maximum overshoot is approximately 25% and the 
maximum settling time is six seconds. The third design is 
basically satisfactory; however, the system design must be 
completed by adding the high frequency poles in L Q . A fourth 

design will now be made to include the high frequency poles. 

. / 

4.8 Addition of High Frequency Poles 

With only minor modification to L Q in the region of one radian 
the third design would be satisfactory from the point of view of 
system response and noise transmission in the low frequency region; 
however, it is far from optimal when noise transmission in the 
high frequency region is considered. The high frequency poles 
necessary to reduce the noise transmission at high frequencies 
will now be added. 

An analysis of the noise transmission problem can be carried 
out in basically the same manner as presented by Horowitz^ for 
linear time-invariant systems. Due to the linearity of the system, 
the system input can be assumed to be zero without affecting the 
results. Referring to Figure 4.1, with the system input zero the 
plant input Z is given by 

Z(jw) = -H(joj) [ N ( jto) + C(jco)] 
which can be written 

1 + H(jw) 

Z(jm) 


Z ( jto) 


= -H (j to) N ( j to) . 
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Rearranging the expression, one can obtain 


Z (juO = 


_ b 0 ( jtlQ 

rjw) 


tea 


Pp(j^ ) 

p eq(j“) 


N (ju) 


+ L 0 (joi) 


(4-10) 


where 


p eq(j“) 


= C(.jai) 
Z(jw) 


and 


L 0 (jw) = P Q (jm)H ( jw) . 


In the low frequency region the characteristics of L 0 are dictated 
by the system response. That is, at low frequencies L 0 must lie 
within acceptable regions in the complex plane up to some frequency 
beyond which | L 0 | may be decreased as quickly as possible. In the 
high frequency region the characteristics of L 0 are dictated by 
Equation (4-10) which describes the noise transmission from the 
noise source to the plant input. To show that it is desirable to 
decrease j L Q j , suppose |l o | were made large in the low frequency 
region. The noise transmitted to the plant input would approxi- 
mately be given by 


Z(jo>)s= 


- 1 
P eq ( 


N(joj) . 


At high frequencies the transmission of the time-varying plant is 
greatly reduced so that |P e q(jo))| becomes small which results in 
large values of noise transmission to the plant input. Thus, |l q 
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should be made as small as possible within the high frequency 
region. 

The reduction of |l q | in the high frequency region is 
accomplished by including high frequency poles in the design of 
L q . A large number of poles cannot be added too low in the 
frequency range for otherwise instability will result. The 
procedure is then to add poles at increasing frequencies as quickly 
as possible while maintaining both a stable system and the proper 
system response. The problem is now to determine at what frequency 

one may stop adding poles . 

Again consider Equation (4-10). At high frequencies |lJ is 
smaller than |P 0 /P e ql so that the noise component of the plant 
input Z can be approximated by the expression 


Z(ja>) 


Ln(» 
Peq(j^) 
PpW 
P eq (i 


N(» 


or 

Z ( j to) L oO w ) N(jw) = H(ja>) N(jw) (4-11) 

PoO) 

From this expression, which is intuitively evident, it is seen 
that poles should be added to L 0 until the desired amplification 
factor between the noise source and plant input is achieved. The 
desired amplification will depend upon the system and the noise 
source. In some cases it may be satisfactory to stop adding 



237 


poles while j H ( j co ) | is larger than unity, while in other cases it 
may be desired to add poles until | H ( j to) | is much- less than unity. 
In this example it will be assumed that poles will be added until 
|H(jw) | becomes less than unity — that is, until | L 0 ( jco) | becomes 
smaller than |P 0 (jo))| as the frequency is increased. 

The approach in making the fourth design is to assume the 
acceptable regions of -L 0 corresponding to the third design are 
also acceptable when high frequency poles are added. Note that 
this does not mean that high frequency poles will merely be added 
to the third design, because the addition of such poles will 
cause an unacceptable decrease in |l o | in the region of five to 
ten radians unless they are placed at extremely high frequencies. 
The high frequency poles must be incorporated as an integral part 
of a new design based upon the assumed acceptable region of -L 0 . 

The polar plot of L Q for the fourth design is shown in Figure 
4.89 together with the acceptable regions and the stability region 
corresponding to the circle criterion. The expression for L 0 (s) 
is given by 


L 0 (s) - 56.16 


s + 5 

s(s + 10) (s + 18) [27 + x ] [fo + X ] [fo + 1 


— + 1 
170 


(4-12) 


The expressions for H(s) and G(s) were not found as a ratio of 
polynominals as was done in the previous designs. Instead, the 
system was simulated by simulating the various components of H and 
G. The expression for H(s) is 
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H(s) + 30.88 gl + . 7s + 10 - - - 
+ 28s + 180 s_ 


+ 1 


+ 1 


4 - 1 


4 - 1 


and the expression for G(s) is 


G(s) = 


1 + L q (s) 


2 . 2s (s + 2 ) 

+ 2.8s + 4 


where L Q (s) is given above. 

In adding the high frequency poles the critical region was 
found to be between five and ten radians per second. As poles are 
added at 20 or 30 radians, lag is introduced between five and ten 
radians which causes -L Q to enter the circular stability region. 
Although entering the stability region does not imply an unstable 
system, system stability is not assured. The magnitude portion of 
the Bode plot of |lJ for the fourth design is shown in Figure 4.90 
with the plot of |L 0 | corresponding to the third design and the 
plot of | P 1 . Note that the third design and the fourth design 
are identical below four radians per second. From four to 25 
radians the magnitude of L 0 for the fourth design is larger than 
the magnitude of L Q for the third design. This increase allows a 
more rapid decrease of |l q | in the region greater than 25 radians. 
The point where |L 0 | for the fourth design crosses |P Q | is at 94 
radians so that poles are not added above 94 radians., 

The fourth design was simulated and the maximum value of |e| 
was determined which is shown in Figure 4.91. It is seen that |e| 
significantly exceeds the specifications in the region of 0.25 to 
1.6 radians. The step response corresponding to a function for 
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Figure 4.90 

Bode Plot of |l o | for Fourth Design, |L 0 |for Third Design and [P 
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which the specification is exceeded is shown in Figure 4.92. 

Note the "dip" in the error as a function of time in the region of 
0.6 seconds. Contrary to what might be expected the "dip" did 
not result in a violation of the specification. The main fre- 
quency component in the "dip" is around 12 radians and the 
specifications were not violated in this high frequency region, 
although an increase in |l d ( was observed in this region over that 
seen in the third design. The violation of specifications in the 
region of 0.25 to 1.6 radians is due to the addition of the high 
frequency poles and is not an unexpected occurrance in a time- 
varying system. 

There are now two possible modifications of |l o | which will 
reduce the value of |e|. Either the magnitude of | L & | may be 
increased in the appropriate low frequency range or the high 
frequency poles may be moved to higher frequencies which results 
in raising [ L 0 | in the high frequency range. The design procedures 
give a formal and systematic approach to increasing |l d | in the 
frequency range where the specifications are violated; however, 
placement of the high frequency poles must be done on a trial and 
error basis . 

The amplification of the noise from the noise source N to 
the plant input is a more severe problem in the higher frequency 
ranges than in the lower frequency ranges. One reason is that the 
noise usually has larger components in the higher frequencies 
than in the lower frequencies. But more important, it is the 
higher frequencies where | P 0 / P e q | is small and |h| is large which 
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Figure 4.92 

Step Response Sample of Fourth Design 
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leads to large amplification of the high frequency noise com- 
ponents as is seen from Equation (4-11). Thus, increasing |l q | 
in the lower frequency range is preferred to relocating the high 
frequency poles to still higher frequencies and thus increasing 
| L 0 | in the high frequency range. 

The fourth design was modified in order to observe the effect 
of both increasing |l o | at low frequencies and relocating the 
high frequency poles. The effect of relocating the high frequency 
poles will first be discussed. The last three poles of |L 0 |, 
Equation (4-12), were moved to 100 radians. This new system will 
be referred to as the first modification of the fourth design or 
simply Mod One. The expression for L 0 (s) corresponding to Mod 
One is 


L 0 (s) 


56.16 


s + 5 


;~(s + 10) (s + 18) 

r S -| 

- + 1 


1 

cn 

+ 

3 


J.1 


L 100 J 



The expressions for H(s) and G(s) is evident from the corresponding 
expressions of the fourth design and will not be given. The 
difference between |l q | for the fourth design and for Mod One is 
most easily observed from Bode plots of the magnitude portion of 
L 0 for the two systems which are shown in Figure 4.93. Although 
the difference in the two designs may not appear to be signifi- 
cantly different from inspection of the Bode plot, the difference 
is actually quite large due to the fact that noise power is 
calculated from a linear integration over frequency. The maximum 
value of | E | for Mod One is shown in Figure 4.94. The specifi- 
cations are violated in the region of .5 to 1.5 radians, but the 






247 


value of | E | is shown in Figure 4.95. Note that the "dip 
which was observed in the error of the fourth design (see Figure 
4.92) is not present in Mod One. This design is unacceptable not 
only because it violates the specifications but also because it 
does not represent the fastest decrease in |l o | that is possible 
in the high frequency region. 

A second modification was made to the fourth design to observe 
the effect of increasing | L 0 j in the low frequency region. The 
modification consisted of increasing the magnitude of |l d | in 
the low frequency range with the high frequency poles unchanged. 
This new system will be referred to as the second modification of 
the fourth design or simply Mod Two. As in the third design, a 
combination of design procedures one and two was used to determine 
the acceptable regions of -L 0 . The shape of the acceptable region 
of -L 0 was assumed to be as shown in Figure 4.89 and the design 
equation, Equation (4-9), was used to determine the proper size 
of the acceptable regions. The acceptable regions of -L 0 for 
the second modification are shown in Figure 4.-96 together with a 
polar plot of -L 0 for Mod Two. The magnitude portion of the Bode 
plots of L 0 for Mod Two together with the fourth design are shown 
in Figure 4.97. Note that in Mod Two the values of |l 0 | do not 
come as close to lying on the boundaries of the acceptable regions 
as does | L 0 [ for the fourth design shown in Figure 4:89. This is 
because it proved to be difficult to have the large magnitude of 
|L 0 | required at 0.628 radians and at the same time decrease |l 0 | 
to the minimum magnitude at 1.25 and 1.88 without violating the 
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Step Response Sample of Fourth Design, Mod 









Figure 4.9 7 

Bode Plot of |L 0 | for Fourth Design 
0 | for Mod 2, Fourth Design and I P 0 
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specifications at 9.42 or even maintaining system stability. Note 
that if the point 0.1818 on the real axis is enclosed by -L 0 , the 
system is unstable when the gain is constant at ten from a time- 
invariant analysis. The plot of -L 0 for Mod Two slightly enters 
the stability circle so that the system stability cannot be 
assured. The system proved, however, to be stable for all time- 
variations tested. 

A slight improvement can be made in fitting -L Q to the 
acceptable regions of -L 0 by moving the pole which is located at 
ten radians to eight radians, moving the zero which is located at 
6.4 radians to four radians, placing an additional zero at 2.3 
radians and placing an additional pole at 1.7 radians. However, 
this small improvement was not felt to be worth the additional 
complexity in the system and will not be discussed in detail. 

The maximum value of the error magnitude for Mod Two is shown 
in Figure 4.98. Figure 4.99 shows a step response for one 
function which corresponds to an extreme value of |e| . Note 
that the "dip" in the error at approximately .5 seconds is larger 
than was observed in the fourth design and is in fact large 
enough to result in a small violation of the specification in the 
region of six to ten radians where |e| reaches a value of 0.6. 

The specifications are also not met in the region of 0.5 to 1.5 
radians . This design would not be considered acceptable for not 
only does it violate the specifications but also stability is not 
assured . 

A fifth design will now be made. It is expected that |l o | 
for the fifth design should lie somewhere between Mod One and 
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Figure 4.98 

Maximum |e| for Fourth Design, Mod 2 
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Figure 4.99 

Step Response of Fourth Design, Mod 2 
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Mod Two, since in Mod One |l q | was too large in the high fre- 
quency regions and not quite large enough in the low frequency 
region, while in Mod Two ) L 0 j was not large enough in the high 
frequency region and appeared to be a little too large in the low 
frequency region. The magnitude of |L q | for the fifth design at 
low frequencies will be chosen to lie approximately half way 
between the magnitude corresponding to Mod One and Mod Two. The 
high frequency poles will then be placed so that a polar plot of 
_L q lies just outside the stability circle in order to assure 
system stability. Note that this approach departs from the design 
procedure but instead draws upon the insight which has been 

gained in the previous designs. 

The expression for L 0 (s) for the fifth design is given by 


L 0 (s) 


75 


s + 6 

7(7 +12.5) fs + 20) fs_ f 

30 





The polar plot of -L 0 is shown in Figure 4.100. For comparison the 
acceptable regions used in the fourth design are also shown in 
Figure 4.100. The relatively large magnitude of L Q in the region 
of nine radians is necessary in order to be able to reduce |l q | 
more rapidly in the higher frequency regions. This can be more 
clearly seen from a Bode plot of |L 0 | which is shown in Figure 
4.101. Also shown in Figure 4.101 are plots of |L 0 | corresponding 
to Mod One and Mod Two. Note that the magnitude of |l 0 | for the 
fifth design lies approximately half way between the magnitudes of 
| l 0 | of Mod One and Mod Two at both high frequencies and low 
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frequencies. The plot of the maximum error magnitudes for the 
fifth design is shown in Figure 4.102. It is seen from Figure 
4.102 that the specifications are satisfied. The choice of L 0 
was fortunate because the value of |e| just meets the specifi- 
cations in the regions of one and 5.5 radians which are the 
critical frequency regions. In the region between two to four 
radians j E | is much smaller than is required by the specification 
which is due to the difficulty in making |l o | as large as required 
in the region of one and six radians, while at the same time 
small enough to just meet the specifications in the region 
between two and five radians . 

The fifth design is very close to the best design possible. 
That is, the design satisfies the specification, but any signifi- 
cant decrease in the magnitude of |lJ in any frequency range will 
either result in a violation of the specifications or an increase 
of the high frequency noise transmission. As has been mentioned, 
the critical frequency ranges are around one and 5.5 radians. If 
the high frequency poles are moved to lower frequencies in order 
to reduce the noise transmission, the specifications will be 
violated at both one and 5.5 radians. Decreasing |lJ in the low 
frequency range will result in a specification violation in the 
region of one radian unless the high frequency poles are moved 
farther out in the frequency range, thus increasing the noise 
transmission. The fifth design is thus selected as the final 
design. 

The expression for H(s) is given by 
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Figure 4.102 

Maximum | E [ for Fifth Design 
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H (s) = 41.25 


(s +2) (s + 6) 


+ 12.5) (s + 20) r-s ,1 rs 1 
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and the expression for G(s) is found from the relation 

G(s) = [l + L 0 (s)l T o (s ) 

Po(s) 


where 


T 0 (s) s(s + 2) 

— — = 2.2 

P Q (s) s 2 + 2 .8s + 4 

A sampling of step responses for the fifth design is shown in 
Figure 4.103 through Figure 4.108. The maximum overshoot was 
found to be 24% and the maximum rise time is six seconds. Observe 
in Figure 4.103 that the "dip" in the error around 0.5 seconds is 
still present but its magnitude has been reduced to the point 
that the specifications are not violated. Note also the high 
frequency variations present in the error shown in Figure 4.104 
and Figure 4.105. These variations are not present in the third 
design, so are thus a result of the decrease in |l q [ in the high 
frequency region. 

The regions of P Q /P e q for the fifth design are shown in 
Figures 4.109 through Figure 4.114. The regions closely resemble 
those found for the third design which are shown in Figures 4.75 
through Figure 4.80. However, the regions for the fifth design 
are somewhat larger than those for the third design, especially 




Figure 4. 103 

Sample Step Response of Fifth Design 
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Sample Step Response of Fifth Design 
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Figure 4.105 

Sample Step Response of Fifth Design 
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Figure 4.106 

Sample Step Response of Fifth Design 
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Figure 4.107 

Sample Step Response of Fifth Design 
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Figure 4.108 

Sample Step Response of Fifth Design 










Figure 4.111 

Region of P 0 /P e q> Fifth Design, to 
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in the critical area around the origin. The addition of the high 
frequency poles in L Q has thus increased the size of the regions 
of P 0 /P e q but tbe basic sba P e is unchanged. 

4.9 Summary 

A system having a time-varying gain which can range between 
values of one and ten has been designed to meet a set of frequency 
domain specifications. The design is close to optimal in the 
sense that noise transmission to the plant input is minimized. 

In the process of carrying out the design a good deal of insight 
about the design procedure and the system was obtained. 

The first comment concerns the specifications. The frequency 
domain specifications which were chosen appear to have been 
reasonable since the resulting system has a well-behaved step 
response that would undoubtedly be considered satisfactory. 
However, the system error as a function of time does not fall 
within the specification envelope shown in Figure 4. 7 which was 
used to arrive at the frequency domain specifications. This is, 
of course, due to the fact which was discussed earlier that error 
functions can satisfy the frequency domain specification and still 
fall outside the time domain error envelope. The time domain 
envelope is exceeded only because of a fast system rise time. 

Since a fast rise time is generally not considered an undesirable 
system property, it appears that the time domain envelope was not 
properly selected in the first place. It appears that the time 
domain envelope should have allowed the system error to deviate 
much more near the origin than was actually allox^ed. If it were 
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desired to have the slower rise time, it would be necessary to 
modify the frequency domain specifications. More study of the 
relationship between the time domain specifications and the 
frequency domain specifications would be necessary, however, 
to determine the specific changes required. 

The approach of using the time- invariant design as a starting 
point appears quite satisfactory. In the example the time- 
invariant design indicated the shape of the regions of -L 0 which 
changed very little during the design. The assumption that small 
changes in L Q cause only small changes in the regions of P 0 /^eq 
proved to be valid for the example. In fact, changes in L 0 that 
were not so small caused only relatively small changes in the 
region of P 0 /P e q as was found in the fourth and fifth design when 
the high frequency poles were added. Although the assumption will 
have to be verified for other types of systems in order to obtain 
any generality, its verification in this example is very promising 
for the use of the regions of P 0 /P e q» 

As was mentioned earlier, an additional advantage to using the 
time- in variant design as a starting point in the procedure is 
that it allows one to readily observe the additional magnitude of 
L q needed to compensate for the time-variations. Figure 4.115 
shows the Bode plot of |l o | for the final design together with the 
plot of |l o | for a time-invariant design which includes the high 
frequency poles. Also shown in Figure 4.115 is the plot of [ P 0 | . 

The magnitude of |l q | for the time invariant design is substantially 
less in both the high frequency and low frequency regions than is 




Figure 4-115 

Bode Plot of |L 0 | for Time-varying and Time-invariant Designs 
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|L | for the time-varying design. However, in the region between 
two and six radians the difference in |l q | between the two designs 
is not large and in fact | L 0 j for the time-varying design could 
probably be reduced somewhat. Thus the time-variations appear to 
have little effect on the system response in the range of two to 
six radians. 

The example also pointed up a serious flaw in design procedure 
one. The flaw, which was found in the second design, is the fact 
that the acceptable regions of -L Q are assumed to be circles 
centered at point one. Such an assumption would undoubtedly lead 
to an overdesigned system with little indication as to just what 
to do to eliminate the overdesign. This is especially true in 
the high frequency region where an overdesign is the most 
undesirable. A problem is also present with procedure two in 
that it is difficult to determine precisely the size of the 
regions of P Q /P eq which means that the precise size of the regions 
of acceptable -L Q is also difficult to determine. It was found 
that a combination of the strong points of the two procedures 
resulted in an improvement in the design procedures. The regions 
°f P 0 /P e q of procedure two was used to determine the shape of the 
acceptable regions of -L Q while |e| of procedure one was used to 
determine how much the regions should be increased or decreased 
for the next design. This combination appears to result in a 
more powerful procedure than either design procedure one or two 
alone . 

It was also found in the example that the addition of the 
high frequency poles in L Q had a large effect on the low frequency 
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error. Thus, it would probably be desirable to include the high 
frequency poles in L 0 from the beginning of the design rather than 
designing L Q at low frequencies before adding the high frequency 
poles as was done in the example. There is, however, an advantage 
in designing at low frequencies first and then adding the high 
frequency poles.' The advantage is that it is possible to examine 
the effect of changing only portions of L 0 . If the poles are 
added early in the design, care must be taken in modifying L Q 
since a change in the high frequency part of L 0 affects the low 
frequency error and vice versa. However, in the example, the low 
frequency portion of L 0 did not affect the high frequency error 
as much as the high frequency portion affected the low frequency 
error. 

It is also of interest to note the type of time-varying gain 
system for which a design is possible. If the time-invariant 
portion of the plant is minimum phase, that is, all poles and 
zeroes are in the left half plane, then it is possible to increase 
the magnitude of L 0 indefinitely without violating the stability 
criterion (see Figure 4.100) so that the specifications will 
eventually be met. On the other hand, if the time-invariant 
portion of the plant is nonminimum phase, then it is not possible 
to increase the magnitude of L Q indefinitely without violating the 
stability criterion. Thus, for such systems it may not be possible 
to increase the magnitude of L 0 sufficiently to satisfy the 
specifications. Therefore, if the time- invariant portion of the 
plant is minimum phase, a design is always possible; otherwise, 
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the possibility of a design depends upon the particular system 
and the specifications. 

Even though a strict application of the design procedure 
itself may not lead directly to an optimal design, the procedure 
does provide a systematic approach to arriving at an acceptable 
design provided' stability can be maintained. In the example the 
insight gained in applying the design procedures allowed the 
design to be carried to a point very close to the optimal. It 
is anticipated that this will generally be the case in applying 
these design procedures. 

The example has demonstrated that suitable frequency domain 
specifications can be obtained and the design procedures can be 
successfully applied to systems with a time-varying gain. 



CHAPTER V 


CONCLUSIONS 

A study of the time-varying system design problem which was 
outlined in Section 1.1 has been undertaken. After a review of 
the previous research in the field of time-varying systems it was 
decided to employ time- invariant compensation and to design for a 
desired transfer function for which procedures are well developed 
rather than attempt to design for an impulse response which would 
require a study in itself. An effort was first made to obtain an 
analytical solution of the system equations from which a design 
procedure could be developed. Unfortunately, the attempt was 
unsuccessful so that efforts were turned to the utilization of 
numerical solutions of the system equations which could be 
obtained from computer simulations. 

It was possible to develop two design procedures which are 
based on the numerical solution of the system equation. The basic 
difference in the two procedures is that procedure one assumes the 
acceptable regions of the nominal loop gain Lq to be circular, 
while in procedure two the shape of the acceptable region is 
determined in the process of the design. Examples were presented 
which illustrated that it was possible to carry out the steps of 
the procedures using computer simulations. 

The procedures as presented do not constitute a complete 
solution of the time-varying system design problem. One 
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difficulty which must be resolved is the determination of 
satisfactory frequency domain specifications. Although it is 
possible to arrive at specifications which intuitively appear 
reasonable, the exact relationship between time-domain and 
frequency domain specifications are unknown. A second difficulty 
is the determination of a general stability criterion for time- 
varying systems. Since necessary conditions for stability are 
known for systems having a time-varying gain, it is usually 
possible to design such systems so that stability is guaranteed. 

(It may not be possible to maintain stability if the time- 
invariant portion of the plant is non-minimum phase; otherwise, 
there will be no difficulty.) Nonetheless, in the design of a 
general system the procedures themselves do not assure stability; 
however, any future stability criterion which is specified in the 
frequency domain can readily be incorporated into the procedures. 

A design example of a system having a time-varying gain was 
carried out in depth. Practical frequency domain specifications 
were developed and a design was obtained which is close to the 
optimal design. At each step the system step responses were 
examined and were compared to the desired step response. It was 
seen that the step responses behaved as one might expect and that 
the step responses of the final design were entirely satisfactory, 
thus showing that the frequency domain specifications were also 
satisfactory. This example demonstrated that it was possible to 
arrive at satisfactory frequency domain specifications and by 
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employing the insight gained in carrying out the procedures, a 
design could be obtained which was close to the optimal design. It 
was seen that in the example the assumption of circular acceptable 
regions of L fl which is made in design procedure one was highly 
inaccurate. In fact, it was found that a combination of procedures 
one and two gave better results than either of the two procedures 
individually. The combined method consists of determining the 
shape of the acceptable regions of L 0 from procedure two, while the 
size of the acceptable regions are determined from procedure one. 

On the whole the example indicated that the design procedures 
could be satisfactorily carried out and would yield acceptable 
results . 

In carrying out the design examples presented in this paper, 
simulations were successfully carried out on analog, digital, and 
hybrid computers. The advantage of the analog computer is its 
ability to solve the system equations quickly; however, it requires 
constant attention by the operator in order to perform the many 
runs required in searching for the maximum E and the region of 
P Q /P e q. Operator action also increases the total run time con- 
siderably. The advantage of the digital computer is that it can 
perform the required runs automatically; however, it requires a 
good deal of time to complete the calculations. The hybrid 
computer is able to utilize the good features of both the analog 
and digital computers. The analog portion of the hybrid computer 
is utilized to obtain a quick solution to the system equations, 
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while the digital portion is used to control the analog and read 
out the results. The hybrid operations are much faster than 
either the digital or analog alone and would be the preferred 
method of performing the calculations if such a computer were 
available . 

An additional area of research is the determination of worst 
case time variations. That is, the determination of those 
variations which result in the maximum value of E and which fall 
on the boundary of the region of P 0 /P e q. At present the method 
of determining the maximum value of E and the region of P Q /P e q 
is to run a large number of possible variations as was done in 
the examples. This approach is successful when there is no more 
than one or two time-varying parameters. However, if the system 
has several time-varying parameters, the number of possible 
combinations of these variations soon become prohibitively large 
to try to determine the worst case variations simply by testing 
representative variations of all possible combinations. If the 
worst case variations were known, there would be no need for a 
search and the time required to carry out the design would be 
greatly reduced. 

It can not be stated the design procedures are ready for 
general application to the design of time-varying systems because 
of the problem of system stability and the problem of determining 
the worst case variations for systems having a number of time- 
varying parameters. However, the procedures can be successfully 
employed on systems having a time-varying gain. In addition, the 
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concepts apply to any time-varying system; so, depending upon the 
case, it may be possible to apply the procedures to more general 
systems. Thus, even though the procedures are limited, they are 
a useful tool in the practical design of time-varying systems. 
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APPENDIX A 


In this appendix it will be shown that if the system remains 
stable, the magnitude of the error can be reduced to as small a 
value as desired by increasing the magnitude of the feedback 
compensation H(jm) . 

Employing the abbreviated notation of Section 3.2, the design 
equation, Equation (3-9), is given by 

P 0 -1 E + HE = -AP" 1 (C 0 + E) (A-l) 

The expressions for E and AP"-*- are 

E = C - C Q 

AP-1 = P _1 - P 0 _1 . 

Substituting these expressions into Equation (A-l) it is possible 
to obtain 


HC + P _1 C = HC 0 + P 0 _1 C o (A-2) 

It is important to understand the meaning of Equation (A-2). A 
block diagram of the system is shown in Figure A.l. 



Figure A. 1 

Block Diagran of System Under Study 
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The system equation is 

GR = P -1 C + HC (A-3) 

Comparing the system equation with Equation (A— 2) it is seen that 
one can write 

HC + P _1 C = GR = HC 0 + P 0 -1 C 0 

That is, both the right hand side and the left hand side of 
Equation (A-2) are equal to the output of the prefilter. This, 
of course, is due to the fact that in Equation (A-2) the plant 
input is considered a constant and the difference in C and C D 
is due only to the difference in the plant P and the nominal plant 
P Q . Thus, the output of the prefilter is always the same 
irregardless of the plant. The system output in turn is such 
that the relation P -1 C + HC is always equal to the prefilter 
output. 

Another point that should be recalled is that in Equation 
(A-3) the plant output is not unique for a given plant input. 

This is due to the fact that the system output is made up of the 
system input and the complementary solution which is uniquely 
determined by the system’s initial conditions .52 in determining 
the step response of a stable system the initial conditions are 
set to zero so that the complementary solution is also zero and 
the output consists only of the forced solution which is uniquely 
determined by the system input. If the system is unstable, 
however, the complementary solution will grow without bound once 
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an initial condition is displaced slightly from zero by noise. 
Thus, if the system is unstable, the complementary solution will 
not be zero so that the system output will not be uniquely 
determined by the system input. Thus, in Equation (A-2) , if the 
system is stable, C is uniquely determined by C Q and vice versa 
since the initial conditions are set equal to zero in the 
determination of the system step response. If the system is 
unstable, however, this uniqueness does not hold. 

In describing the system it is convenient to borrow notation 
from classical linear time-invariant system analysis. Let the 
system type number of a time-varying system be the same as the 
type number if the plant parameters were time- invariant A 1 That 
is, if the plant parameters were time-invariant, then the plant 
of a type 0 system has no poles at the origin, while the plant 
of a type 1 system has a pole of order one at the origin, the 
plant of a type 2 system has a pole of order two at the origin, 
and so forth. 

Consider next the relationship between the nominal plant and 
the desired step response. The desired step response C D must be 
chosen such that it can be delivered by the nominal plant P Q as a 
result of a reasonable plant input. That is, it must not be 
necessary that the plant input contain impulse functions or 
derivatives of impulse functions in order for the nominal plant 
to deliver the desired step response when the plant is the nominal 
plant. This will only be true if the input to the plant is 
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physically possible. Thus, if the system is stable, the input to 
the nominal plant P Q must be bounded. In addition, if the system 
is type 0, e(t) will settle to a constant and if the system is 
type 1 or higher, e(t) will settle to zero. Therefore, for type 
0 systems the Fourier transform of the plant input will only have 
a pole at frequency zero and no other poles, (i.e. , all the poles 
of the Laplace transform of z(t) will be in the left half plane 
except for one pole at the origin). For type 1 systems or higher 
the Fourier transform of the plant input will have no poles (i.e., 
all poles of the Laplace transform of z(t) will lie in the left 
half plane). The Fourier transform of the input to the nominal 
plant is given by 



which is independent of the feedback H; thus, the magnitude of 
p ~lc is independent of H and is bounded for all co except in 
type 0 systems where T 0 ~'*'C 0 has a pole at w = 0 . 

The desired objective is to show that the magnitude of the 
error goes to zero as the magnitude of H goes to infinity provided 
the system remains stable. Since the system is assumed stable and 
for step responses the system initial conditions are zero, Equation 
(A-2) uniquely describes C in terms of C 0 and vice versa. Equa- 
tion (A-2) can be written in the form 



Returning to the full notation, one obtains 
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E(jw) = C(jco) - C 0 (joj) = 

00 

Cp(juj)Po 1 ( joj) - J P“1 ( jcj, j y) £ ( 3Y ) dy. (A-4) 

H(ju) ““ H(ju) 


Note that the nominal plant is time-invariant and can thus be 
represented by its transfer function which is denoted P 0 - ^(jm). 

From the previous discussion it is seen that the first term on the 
right hand side of Equation (A-4) goes to zero as |h| goes to 
infinity except for type 0 systems at id = 0. This is because 
P 0 _1 C 0 is independent of H and is bounded except for the one case. 

It is not so obvious, however, that the integral term goes 
to zero as | H | goes to infinity since the system output C is 
dependent upon H. It will be shown by contradiction that the 
integral term does in fact go to zero as the magnitude of H 
increases without bound. Suppose that the magnitude of the 
integral term does not go to zero as |h| goes to infinity. That is 


lim 
] H |-*» 


/ p_ 1 (j“,jy) dy 


H(jco) 


> K. 


CA-5) 


Observe that H(jm) is not a function of the variable of integra- 
tion and thus just represents a complex number that grows in 
magnitude with the limiting process. Interchanging the order of 
the limiting process and the integration procedure and noting that 
P - ^ ( jw > jy) is independent of the limiting process since it is 
independent of H, one can write 
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00 

/ P -1 (jw,jY) 

—CO 


lim C(.jy) dy 
| H 1 H(jtu) 


> K. 


(A- 6) 


In this form it is evident that the magnitude of C must increase 
at least as fast as the magnitude of H, for otherwise, the 
magnitude of C/H would tend to zero and the magnitude of the 
integral would also tend to zero. Thus, C must have the property 

lim j c | A ] H | b , b > 1, A > 0. 

|h|-~» 


Observe, therefore, that 

lim Ip _ r 

IhI-w |c 1 


AH 


Thus, from Equation (A-4) , one now obtains 


lim 


C - C Q | | H I = 


uu 

/ p-h^.v) 


or 

lim A | H | b | H | < /]p- 1 (jo ) ,jy)| A l H | b d T- 

| H I"* 00 —oo I H I 

Dividing by |H| b+1 , the expression becomes 

00 

4 s L |F ' 1(3 “’ 3T)I iHh»m dY ' ' 

Since the term A/|h| goes to zero as the magnitude of H goes to 
infinity, the integral is zero or 


A < 0 



which is a contradiction. The original supposition that the 
integral term of Equation (A-4) does not go to zero is thus false, 
and it can only be concluded that the integral term does indeed 
go to zero as the magnitude of H goes to infinity. 

Since both terms on the right hand side of Equation (A-4) 
go to zero as the magnitude of H increases, the magnitude of the 
error must also go to zero. 

In summary, it has been shown that as the magnitude of the 
feedback compensation H is increased without bound, the magnitude 
of the error tends to zero provided system stability can be 


maintained. 



APPENDIX B 

DERIVATION OF I (jy) 

In this appendix the expression for I (jy) of design example 
two is derived. The expression for I Q (jy) is given by 

co 00 i 

I (jy)= J‘ p" 1 (jy,jtu)c o (j^)dtu - J' P^ (jy,ju))du) 

° -CD -oo 

which is equivalent to the Fourier transform of the plant input 

when the plant output is equal to c^(t) minus the Fourier transform 

of the nominal plant input when the output is equal to c^(t). 

Referring to design example two it is seen that if the plant output 

is c (t) then the plant input is given by 
o 

y(t) = g(t)c' o (t) + f(t)c Q (t) , 

and if the input to the nominal plant is C Q (t) its output is 

y(t) = 0.505c (t) + 0.175c (t) . 

J o o 

Thus the expression for I (jy) can be written 

o 

I (jY)= J {(g(t)-0.505)c* o (t) + (f (t)-0.175)c o (t)j e jYt dt . 

° -oo (B — 1 ) 

Substituting the expressions for g(t), f(t) and C Q (t) from design 
example two into this equation one obtains 

I o (j y)= J 1^-0.505 4- G 2 (l-e" b(t " T 2 ) )u(t-T 2 )] 

[4e _1 ' 4t cos(1.414t) - 3 . 96e _1 ,4t S in (1 .414t) ] u(t) + 

[f i - 0.175 + F 2 (l-e' a( ' t " T l ) )u(t-T 1 )] 

^2.83e _1 " 4t: sin(1.4.4t)u(t)j j- e' jYt dt . 

Note that the evaluation of this integral depends upon the values 
of and t 2 . Let 

G 3 = G x - 0.505 
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and 

F 3 = F x - 0.175 . 

In order to simplify the expression the following functions are 
defined : 


A^jy)- £ G 3 [4e' 1 - 4t cos(1.414t) - 3 . 96e _1 - 4t sin (1 .414t) ] e‘ jYt dt 

A 2 ( j Y) = r G 2 (l - e” b< ' t ’ T 2' ) ) [ 4e 1 ,4t cos (1 .414t) - 

0 -1 4t 1 -iyt 

3.96e ’ sin(1.414t) J e dt 

A 3 (jY)= J G 2 (l - e" b< ' t ' T 2' > ) [ 4e 1 * 4t cos (1 .414t) - 

T 2 3.96e‘ 1 ' 4t sin(1.414t) ] e‘ jYt dt 

CO 

B 1 (jy)= F 3 2.83e 1 ' 4t sin (1 .414t) e - )Yt dt 

B 2 (j Y )= J F 2 (l - e" a(t ‘ T l } ) 2.83e" 1,4t sin(1.414t) e" jYt dt 

B (j Y )= J F 2 (l - e -a(t-T i ) ) 2.83e' 1-4t sin(1.414t) e‘ jYt dt 
T 1 

Evaluating these integrals one obtains: 


A 1 (jY)= G 3 


A 2 (jY)= G 2 


4jy 

(jy + i • 4) 2 + 2 

r 4j y _ 4 e T 2 (jy + b) 

(jy + 1.4) 2 + 2 (jy + 1.4 + b) 2 + 2 


-(jy+1.4)T f 4 j ycos (1 .414T2) -3 . 96 ( jy+2 . 8) sin (1 .414T2) 

A 3 (jy)= G 2 e 2 -2 

I (jy + i.4) + 2 

4 ( j yfb )cos( 1 . 414 T 2 )- 3 . 96 (jyf 2 . 8 +b )sin( 1 . 414 T 2 ) 

(jy + 1.4 + b) 2 + 2 


B 1 (jy)= f 3 


(jy +1.4) +2 

r 4 


aTi 

4 e 1 


2 " 2 
(jy +1.4) +2 (jy + 1.4 + a) +2 


B 2 (jy)= f 2 
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B.J ( j "V) = F 0 e 


( jyfl .4) f 4cos (1 .414t^)+ 2 . 8 ( j Y+l .4) sin (1 . 4t^) 

l ( j y + 1.4) 2 + 2 

4cos (1 . 414t^)+2 . 8 ( j y+1 • 4+a) sin (1 . 414t^) 

(jy+1.4+a) 2 +2 J 


I (j y) is thus given by the following expressions. 
For t < 0 , t 2 <■ 0 

I q ( j Y) = ( j Y) + A 2 (jy) + B^jy) + b 2 ( j y) 

For t < 0 , t 2 > 0 

I ( j Y) = A 1 (jy) + A 2 (jy) + Bj ( j y) + B 3 (jY) 

For > 0 , t 2 < 0 

I ( j y) = A 1 (jY) + a 3 ( j y) + B^jY) + B 2 (jy) 

For t > 0 , t 2 > 0 

1 ( j y) = A^jY) + A 3 (jy) + B x (jY) + » 3 (jY) 


Q 



APPENDIX C ' •' 

DERIVATION OF P (juo)/P (j a ,) 
o eq 

In this appendix the expression for P (j(ju)/p (juo) of design 

o eq 

example three is determined. Referring to Equation (3-55) 

P (j(Ju)/P ( j < jo) is seen to be given by 
o eq 

P o (i-) P (j«0 00 l 

p (jm) " FIFO S P ' (j».jY)C 0 «^ d Y ■ 

eq o -oo 


The transfer function of the nominal plant is given by 

1 


P (j(t) = 


0.505 ( j oo) + 0 . 175 j oo 


and C Q (juj) is given by 


C (j«0 = 


joo J^(ju)) 2 + 2.8jw> + 4 J 


Let 


Y(joo) = J P' 1 (jou,jY)C o (j Y )dY . 

-00 

This integral is the Fourier transform of the plant input when the 
plant output is c^(t) . Thus the integral can be expressed as 

CO 

Y ( j to) = J / g(t)c' o (t) + f(t)c Q (t)l dt . ( C— 1 ) 

Comparing Equation C-l with Equation B-l of Appendix B it is seen 
that the two are the same with the exception of two constant terms 
in Equation B-l. Thus the development given in Appendix B applies 
directly to obtaining Y(juo) with the exception that is set equal 
to and F^ is set equal to F . That is 

Y o (j«0 = I c> (j UJ ) 

where 


G 


1 



and 
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The expressions for P (joi)/P (jtr) are 
r ’ o eq 


For t < 0 , t < 0 


P 0 (jw) P 0 (j«3) 
= C 0 (jUj) 


A x ( j oo) + A 2 ( jto) + B 1 (jiu) + B 2 ( j cu) J 


For T < 0 , t 2 > 0 


P Q (j«)) P Q (jti)) 
= C o (jo)> 


[A 1 (jua) + A 2 (jio) + B 1 (ja)) + B 3 (jU))J 


For t > 0 , t 2 < 0 

P (ju>) P q ( j tt>) 

p eq (M = MM 


A 1 (jd)) + A 3 (jo)) + B x ( j cu) + B 2 ( j cu) J 


For T > 0 , t 2 > 0 

P (ju>) P (jcu) 

o _ _o 

P eq (j^) C o (jcu) 


[ A l (jo>) + A 3 ( j cu) + B q ( j cu) + B 3 (juj)] 


where the functions A^ju)), A 2 ( j cu) , A 3 ( j cu) , B^jcu), B 2 ( j cu) , and 
B ^ ( j cu) are given in Appendix B with G q equal to and F^ equal 
to F„ . 



APPENDIX D 


This appendix discusses the computer simulation of the design 
example presented in Chapter 4. 

The simulation was carried out on an IBM 1130 digital computer 
using the Continuous System Modeling Program or CSMP for the 1130. 
The computer had a core storage of 16K and a memory access time of 
four microseconds. The CSMP program simulates an analog computer 
with the input in the form of functional blocks such as an integra- 
tor, summer, multiplier, etc. The system to be simulated must 
thus be represented by a block diagram in much the same manner as 
is done in setting up an analog computer simulation. The program 
is limited to a maximum of 25 integrators and a total of 85 
functional blocks. This size proved to be adequate for the systems 
simulated in the example. The method of performing the inte- 
gration in the program is by a second order Range-Kutta which is 
carried out at one half the integration interval specified by the 
user . 

Modifications to the CSMP program were necessary in order to 
save and store the plant input and output signals as well as to 
calculate the necessary Fourier transforms. The algorithm given 
in Reference (48) was used in the subroutine for the Fourier 
transform computation. Modifications were also made in the program 
for selecting a different time-varying function at each run. The 
function parameters could either be selected at random from a . 
random function generator or, if desired, could be chosen by the 
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operator and read in from cards. Provisions were made for 
calculating and printing the error E(jto) and the ratio 
Pod 00 ) / p eq(j w ) at se l ecte d values of to. In addition, the actual 
signals themselves could be printed out if desired. As an aid in 
determining the regions of P Q /P e q, a plotting scheme was devised 
whereby the points of P Q /P eq were stored for each calculation 
and after all time-variations had been run, the points were 
"plotted" by the printer for each value of to. This technique 
eliminates the necessity of manually plotting each value. 

The calculation of the Fourier transform of the various 
signals requires special consideration. In the definition of the 
Fourier transform, the limits of integration are from plus to 
minus infinity. Obviously, the numerical integration cannot be 
performed over such a range; however, if the function to be 
integrated is zero outside a finite range, the integration need 
only be performed over this finite range and the numerical 
integration can be performed satisfactorily. Consider the error 
function 

e(t) = c(t) - c Q (t) . 

At t less than or equal to zero, c and c 0 are both zero so that e 
is also zero. If there is an integrator in the plant, the final 
value of c will equal the final value of c 0 , so that the final 
value of e is also zero. Practically, one finds that e is 
essentially zero after some finite time T. Thus, the numerical 
calculation of E(jto) is accomplished by integrating over the 
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interval zero to T. Difficulty is encountered in the numerical 
evaluation of the Fourier transform of the step response c(t). 

Note that c is zero for t less than or equal to zero but goes to 
some finite value C m as t goes to infinity. Thus, the numerical 
evaluation of C(jm) cannot be made from c(t) directly. This 
difficulty can be circumvented by defining a new function c(t) as 

c(t) = u(t) - c(t) 

where u(t) is the unit step function. After some time T, c(u) 
is essentially equal to c so that c(t) is zero for t greater than 
T. Also, c(t) is zero for t less than zero so that the numerical 
Fourier transform of c(t) can easily be calculated by integrating 
over the interval zero to T. The Fourier transform of u(t) for 
cu greater than zero is known to be 

u(t) = i- 

so that the Fourier transform of C can be obtained from the 
relation 

C(jm) = Cjf'j^u(t)j -<f[c(t)] 

= — - C( joi) 

The function c(t) can easily be obtained in the simulation from 
which C(jw) and C(ju)) are then calculated. The calculation of 
the Fourier transform of the plant input presents no particular 
problem since in the example the plant input signal goes to zero 


as time goes to infinity. 
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The integration step size in the program must be selected 
by the operator and is unchanged throughout the run. A number of 
step sizes were examined in order to determine the proper values 

t 

which would yield as rapid a solution as possible and still 
maintain accuracy. The size selected was 0.01 seconds. Tests 
were run on time-invariant systems to check the accuracy of the 
Fourier transform calculation and it was found that the computer 
computations error was well within 5% up to 30 radians. As an 
additional test, each time-varying design was run with the gam 
held constant at values of 1, 1.818, and 10, and the results 
compared against analytical calculations. This procedure not 
only verified the accuracy of the calculation but also insured 

against an error in the simulation. 

Despite efforts to reduce computation time, the calculations 
were fairly lengthy. It required from seven to ten minutes of 
computer time to calculate ten seconds of simulated system time 
and an additional two minutes were required to compute |e| and 
Po/p e q for the frequencies of interest. A complete design 
evaluation consisting of from 250 to 300 runs required from 35 


to 40 hours of computer time. 



